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INTRODUCTION 


Literature and data sources 

The literature on the properties and spectra of ocean waves and 
water waves that have been generated by the wind Is substantial. Beginning 
with "Ocean Wave Spectra" and continuing through Kinsman (1965), Neumann 
and Pierson (1966), Phillips (1966), and Hess, Hidy and Plate (1969), various 
texts and review papers have presented summaries and interpretations of 
the data and their analysis. Additional material subsequent to these analyses 
has become available from observations in wind tunnels made by Toba* , 
Sutherland (1967), and at New York University, from measurements in the 
Black Sea by Leykin and Rosenberg (1970, 1971) and from numerous investi- 
gations of the gravity wave spectrum for open ocean conditions. 

Most of this work has concentrated on measuring the time vari- 
ation of elevation at a point, but other studies such as Cote et al.(1960), 
Schooley,** Schule, Simpson and DeLeonibus (1971), and Dobson (1970) 
have reported on measurements of elevation as a function of position either 
on a rectangular grid or along a line. Cox (1958) has measured high fre- 
quency slopes in the upwind -downwind direction for very short fetches as 
a function of time. 

For an appropriate spectral definition and a modest range of speeds, 

8 3 

the wave number spectrum can range from 10 (or so) cm at a peak to 
10 - ^ (or so) cm^ near the viscous cutoff region and the wave numbers can 

^Personal communication of data. 

**Personal communication of data. 
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range from 2 tt 10 cm (a wavelength of 1 kilometer) to 2 tt 10 cm 

(a wavelength of 0.1 cm). There are thus more than thirteen orders of 

magnitude in spectral values to be covered in attempting to understand 

the spectrum of a wind- roughened sea which in turn covers about six orders 

of magnitude in the wave numbers involved. 

Two reports by Martin (1966, 1969) were particularly helpful at 
the start of this investigation, Martin attempted to recover wave number 
spectra from the limited amount of data available to him so as to apply 
them to problems in underwater sound. These efforts resulted in an indi- 
cation of strong variation with wind speed in the capillary wave portion of the 
spectra and an indication of the need for some transitional region to con- 
nect the gravity wave portion to the capillary portion. 

Measurements under natural cond itions 

This entire set of published material is not totally consistent for 
many reasons. Wave recording systems on ships and on offshore towers 
such as those once at Argus Island have to be rugged enough to withstand 
stormy seas and therefore cannot respond to, and measure, the high frequency 
part of the spectrum. The Tucker shipborne wave recorder (Tucker, 1956) 
more or less cuts out about « = 2tt/4 and various artifacts have had to be 
employed to handle the calibration and digitization noise problem that de- 
velops when wave records taken by this instrument are analyzed as in 
Moskowitz, Pierson and Mehr (1962, 1963). The laser altimeter data of 
Schule, Simpson and DeLeonibus (1971), for example, cuts out a wave 
number of 0.36 m 1 corresponding to a wavelength of about 17 meters. 
Additional open ocean wave data given by Saenger (1972) provide data of 
value on the frequency spectrum from f = 0.35 to 0.7 hertz by means of 
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data obtained with a very accurate version of the NIO heave, pitch and 

! 

roll buoy. 

Although many open ocean, marginal sea, lake, reservoir, and 
wind-water tunnel records of waves generated by the wind have been 
spectrally analyzed and reported in the literature, not too much has been 
said about the calibration and frequency response of the recorder systems 
that were used. A notable exception is the field study by Kinsman (I960) 
in which every step of the calibration and recording process was checked 
and described. These data were of such high quality that they were used 
to document the non-Gauss ian probability density function for wave elevation. 
These non-Gaussian properties are important enough to consider in refined 
theories of wave properties, but it is still important to determine the wave 
spectrum as the resolution of the variance of the wavy surface into fre- 
quency and direction bands, or into vector wave number space. 

Measurements in wind-water tunnels 

Measurements of the high frequency, high wave number, portion 
of the wave spectrum have been carried out in wind-water tunnels under 
more easily controlled conditions and with more sensitive instrumentation. 

It is, of course, then necessary to assume that these short waves as 
measured in a wind tunnel superimpose on the long waves measured on 
open water. The formulas for the generation of wave spectra due to Miles 
(1957, 1959a, 1959b, 1962) and Phillips (1957, 1966) involve the spectra 
of the advected turbulent wind fluctuations. These spectra might differ from 
one facility to another and from what occurs on the open ocean, where, for 
example, there are scales of gustiness that are long compared to the time 
constants involved in the generation and decay of capillary waves. This 
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does not seem to be what happens when differences in the nature of the 
winds over the ocean compared to winds in wind-water tunnels are ac- 
counted for in a reasonable way. 

The gap in the data 

The study of these various sources of data, along with recognition 
of the problem of instrument calibration for some of the data, led to the 
realization that a portion of the frequency band (or wave number band) 
that had to be covered was not well measured in either the field data or in 
the laboratory data. The laboratory data do not yield low enough fre- 
quencies because of the limited fetches, and the field data do not extend to 
the higher frequencies. 

This data gap has, however, just recently been removed by the 
work of Leykin and Rosenberg (1970), who made measurements in the Black 
Sea covering frequencies from 1 to 7 hertz for a variety of wind speeds 
with a recording system whose calibration was known and checked in the 
laboratory. The recording system was such that the high frequencies were 
pre-whitened, which resulted in more reliable results. 

These measurements, plus wind tunnel measurements, make it 
possible for the first time to obtain theoretical spectra that are consistent 
with the results obtained in the classical study of sea surface slopes as a 
function of wind speed made near the Hawaiian Islands by Cox and Munk 
(1954). The slope statistics of Cox and Munk have been verified over a 
limited range of wind speeds from spacecraft photographs of glitter patterns, 
the difficulty being that such patterns are usually available only in the trade 
wind areas of the subtropics for a modest range of wind speeds. 
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An Idealized spectrum for idealized conditions 

The concept of a wind generated sea whose properties are solely 
determined in a neutrally stratified atmosphere by the direction and speed 
of the wind that has blown a sufficiently long time over a sufficiently large 
area of the ocean has been questioned. Apart from the fact that large 
enough areas and long enough durations do not occur for the high winds, 
and that the waves are modified by swell and left-over dead seas for low 
winds, it is believed by some scientists that the spectrum of the waves will 
continuously grow toward low frequencies if a constant wind continues to 
blow over a long fetch (Phillips, 1966). 

The main reasons for questioning this belief are the characteristics 
of the trade wind seas, which seem to be in equilibrium with the winds over 
vast areas, and the fact that the lowest frequencies observed in ocean 
wave spectra correspond to periods of about 24 seconds and never seem 
to get lower. Such a continuous growth toward lower frequencies ought to 
produce these low frequencies, if they were possible. 

For these reasons, in this study, a fully developed wind wave spec- 
trum that is solely a function of u^ and U 19 as related to. u^ through 
the theory of a logarithmic wind profile in a neutral atmosphere, is given. 

The actual spectrum for a particular area of the ocean will depart from this 
theoretical spectrum for a variety of reasons. 

For the lower frequencies (and longer wavelengths), the spectral 
components require a relatively long time to be generated and lag the build 
up of a local wind. Once they are generated, they persist as a M dead n sea 
and propagate into areas of lighter winds as swell. Moreover, there is in- 
creasing observational support from both the JONSWAP study (Barnett, 1970) 
and from spectra obtained by the Tucker shipborne wave recorder, Hoffman* 


^Personal communication. 
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(unpublished), that nonlinear interactions as proposed by Phillips (I960, 

1966) and studied by Hasselmann (1962, 1963) and Mitsuyas u (1 968) 
modify the spectra by exporting spectral energy from the central fre- 
quency bands both to low frequencies so that the spectral peak becomes 
more pronounced (and perhaps shifts to lower frequencies) and to high fre- 
quencies where spectral saturation causes breaking and the destruction 
of the organized wave energy. 

Our impression at present is that the fully developed spectrum 
given by Pierson and Moskowitz (1964) probably applies as long as the wind 
continues to act on the waves because the middle band will continuously be re 
generated. However, once the wind decreases, the spectra will die down 
in the middle gravity band and increase in the low frequency band. The 
spectral energy transported to the high frequency gravity band will be 
destroyed by breaking. 

Measurements in the high frequency gravity band and on into the 
capillary band made in a variety of ways with a variety of instruments for 
both field conditions and in wind-water tunnels all show that the spectrum 
at high frequencies is strongly variable as a function of wind speed. There 
is no equilibrium range in the capillary spectrum. There is a range of 
frequencies in the gravity wave range that have an equilibrium value, but 
this range is wind-speed dependent and shifts toward lower frequencies 
with increasing wind. 

For a given area of the ocean at a given time, the spectrum of the 
waves in that area will differ from the spectrum given here in a number of 
ways. The low frequency, low wave number, portion will often correspond 
to higher waves with swell bands when the winds are very light. If the winds 
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have died down from, say, 20 m/s to 10 m/s over a large area, the 
spectrum will develop, with time, a sharper low frequency peak and a 
deficit at middle frequencies. The higher frequency gravity portion will 
be in equilibrium with an equilibrium range over a portion of the band 
and with the wind over another portion, probably with the wind correspond- 
ing to the properly averaged synoptic wind. 

The capillary portion of the spectrum will be in practically instant- 
aneous equilibrium with the gusts and lulls of the turbulent wind over 
areas with dimensions of 10’s of meters. These scales of wave response 
were described by Pierson and Ivloore (1972) in interpreting the variability 
of radar backscatter measurements made under the MSC earth resources 
program. The spectrum in the capillary range will also have angular 
properties that respond to the lateral variation of the turbulent wind about 
the average wind direction. 

It is thus possible to imagine that the actual spectra of wind- 
roughened seas will have shapes roughly like the ones proposed at low 
wave numbers but depart from the theoretical form because of the variations 
in time and space of synoptic scale winds and because of swell. At middle 
wave numbers, the spectra will be in equilibrium with the local wind as 
averaged over 20 minutes or so. At high wave numbers, the spectra will 
be in practically instantaneous equilibrium with the turbulent gusts and 
lulls of the wind. 

Analysis procedures and the fundamental hypothesis of this investigation 
The analysis procedures used in this study are the standard pro- 
cedures for the analysis of stationary time series data in terms of spectra, 
where the statistical term "time series’ 1 also includes the extension to space 
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series and two-dimensional surfaces. Fast Fourier Transforms were 
used, if needed, to sharpen up the analysis. When appropriate, the analysis 
procedures are given in greater detail in the various sections of this study. 

The fundamental hypothesis in this study is, however, very differ- 
ent from the hypotheses made by many other scientists in the study of 
similar data. The hypothesis is simply that once the effect of fetch, as 
some multiple of the wavelength involved, is overcome so that the spectral 
component is generated, its value is solely determined by the speed of the 
wind. This hypothesis assumes, then, that waves with lengths of 10 meters 
do not change properties for fetches of one kilometer or 100 kilometers. 
Waves with lengths of 10 cm do not change properties for fetches of 100 
meters or 10 kilometers. This hypothesis, as tested against a wide 
variety of data sources, makes it possible to demonstrate that the higher 
wave numbers, and higher frequencies, are strongly dependent on wind 
speed. 

Analysis procedures that normalize the frequencies in the spectra 

to f/f max where f max is the frequency of the spectral peak or that use non- 

dimensional frequencies of the form u f/g and that then use non-dimers ional 

2 

fetches of the form gF/u^ relocate the same true frequency in different 
places in the normalized spectra and in the analysis. Such procedures may 
lead to strange conclusions such as the one that the spectral constant in 
the so-called equilibrium range decreases with increasing fetch. 

Wave spectra and radar backscatter 

The reason for this study is that there is a need for a better 
definition of the spectrum of a wind- roughened sea surface for use in the 
interpretation of radar sea return measurements. There is also need 
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for Improved theories of radar sea return, especially for measurements 
near vertical incidence and out to 30°, or so, off the vertical. 

Wright (1968), for example, used the available literature to 
reason that both the high frequency gravity waves and the capillary waves 
were saturated (with perhaps different spectral constants) and did not 
change with wind speed, and concluded from Bragg scattering theory that 
the measurements of radar backscatter should also saturate at a rela- 
tively low wind speed. As will be shown, the high wave number region 
of the wave spectrum varies strongly with wind speed and Bragg scattering 
theory can thus be used to relate these wind speed dependent wave number 
regions to radar backscatter measurements. Radar sea return measure- 
ments should not saturate with increasing wind speed, as far as the 
spectrum is concerned. Whitecapping and foam are not treated in this 
analysis and could alter the real situation at high wind speeds. 

Beginning with the work of Chia (1968), it has become clear that 
a better definition of the wave spectrum is required for improved theories 
of sea return. The need for an improved wave spectrum becomes parti- 
cularly acute in the recent theoretical work of Jackson (1973), which ex- 
tended the work of Chia. This theory for radar sea return requires a 
vector wave number spectrum with fourth moments in both t and In (to be 
defined later), and the spectrum proposed by this study has these moments. 
A final task in this Langley AAFE Program as described by Pierson, 
Jackson, Stacy and Mehr (1972) is to evaluate the equations given by 
Jackson (1972) using the spectra given in this report so as to determine 
theoretical sea return values for the near vertical. 
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DEFINITIONS AND NOTATION 


Spectral forms 

Given a large area of the ocean in which the sea surface ele- 
vatlon is known as a function of x* and y* as q(x*, y*) then the covariance 
function can be defined as equation (1. 1), as L becomes large. 

x*+L/2 y*+L/2 

R ( x > y) = — ^ J ^ t|(x* f y*) n(x*+x,y* + y) dx*dy* (1.1) 

k x*— L/2 y* -L/2 

It is assumed that the wind is blowing in the positive x-direction. 

The Fourier transform of R(x, y) defines the wave number 
spectrum S or S(k) where |k| = (£ +m f 

oo oo 

S J ^ R(x, y) CO s { Ax + my) dx dy (1.2) 

F - oO O 

under the assumption that S*(j t,m) = 0 if i < 0 because of the assumed 
wind direction. 

With l = k cos $ and m = k sin$ , it follows that S*U, m) d4dm 
can be transformed to 

S*(kcos$, k s in $) k dk d$> = S *(k, $) k dk d$> (1.3) 

= S (k, $) dk d<5 

Furthermore, if positive x is in the wind direction and if 
S(£,m) = 0, l < 0, 

tt /2 

/ S(k,*)d* = S(k) (1.4) 

-tt/2 
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The Integral from 0 to oo of S(k) is the variance of the random process 
and equals R(0,0). Equations for S(k,$) will be given below according 
to the above definitions. 

Now, if the sea surface can be represented by a stationary 
linear Gaussian process, it is possible to discuss the spectra of slopes 

and curvatures. The spectrum of 8r|(x, y)/8x, for example, is given by 

2 2 
1 S(Z, m) and the spectrum of 8t](x, y)/8y is given by m S(i,m) . The 

spectrum of the total slope of the sea surface is thus 

S {i,m) = U 2 + m Z ) S(i,m) (1.5) 

s 

which transforms to 

S fl (k) = k 2 S(k) (1.6) 

upon integrating out and for example 

(T 2 = o- 2 + cr 2 = cr 2 + o- 2 = f k 2 S(k)dk (1.7) 

x y u c J 
3 o 

should be related to the measurements of slope variance made by Cox 
and Munk (1954). 

It can be shown that a similar measure of curvature is 

given by 

S c (k) = k 4 S(k) (1.8) 

This function will also be discussed along with the validity of the 
assumption that the sea surface can be thought of as a linear Gaussian 
process for such operations. 
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Wave number - frequency relationships 


The relationship between frequency and wave number for 
the full range of frequency and wave number is given by equation (1.9). 


or by 


whe re 


and 


U = w(k) = (gk + 1 k 3 ) 2 


^ 2 


CO = GO 

m 2 


k k 3 n2 
k + , 3 


L m k 


m 


k m = 


w m = "^8 <gP /T ) 


and where t is the surface tension. 

The inverse of equation (1.10) is equation (1.13). 


(1. 9) 


(MO) 


( 1 . 11 ) 


(M2) 


k ~ k(oo) = 3 co /gF (go /to ) 


(M3) 


1 1 

where F(<o/u> ) = 1 + [1 + B(w) + (B(tu)(2 + B(oo)))^]^ 

X ^ 

+ [1 + B(w) - (B(«)(2 + BM))2 ] 3 (1.14) 

4 

and where B(oo) = 54 ( 00 / 00 ^) as given, for example, in Neumann and 
Pierson (1 966). 

Some questions have been raised about the validity of using 
these relationships in the presence of possible currents in the water. 

A time history, r)(t), yields a spectrum S(oo) that it would be desirable 
to transform to S(k) by means of 
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S(k) dk - S(w(k))^|M dk 


(1.15) 


If the waves are being translated on a true current in the wind 
direction, the frequency at a fixed point would be higher than the true 
frequency that goes with a given wave number* However, there can be 
currents in the water due to Stokes' drift as shown by Chang (196^) for 
gravity waves and yet the relationship between frequency and wave 
number would still be the one given above. Equations (1.9) and (1.13) 
will be used to transform spectra back and forth as needed from wave 
number space to frequency space. The various data sources support 
the validity of this procedure, as will be shown. 

It should be pointed out that systems of measuring and analyzing 
wave elevation along lines, or on a rectangular grid, do not yield spectra 
that are easily related to S(k). For example, the measurement of sea 
surface elevation along a line parallel to the wind direction yields a 
spectrum of the form 

oo 

SU*) = / S{S,,m)dm (1.16) 

- OO 

This expression is not easily related to S(k), and various complex ex- 
pressions result when S ( j&, /7Z ) is wind speed dependent and defined in 
different ways over various ranges of k. Schule, Simpson and 
DeLeonibus (1971) have studied how spectra of this form in the gravity 
wave range can be related to directional spectra under various assumptions 
for the angular spread of the spectrum. 
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AN ANALYTIC REPRESENTATION OF THE SPECTRUM 
OF A WIND ROUGHENED SEA 

For a friction velocity, u ; ,_ greater than some minimum value, 

= 12 cm/seci the spectrum can be defined by the following equations. 

S(k, ft) = S{k) F(k,ft) ; 0<k<oo, - f < $ < £ 

S(k, ft) = 0 ; 0 < k < oo . ^ < ft < it , -tt < ft < - ^ 

F(k, ft) = -J- (cos ft) [1 - e" g /2k ] 

JTT 

+ I[1 - I + A(cosft) 2 ]e- g2/2k2 ( U(u ^ )4 

TT ^ 

where 0 < A < Z and A could be a function of k and u . if needed. 


Sj(k) ; 

0 < k < k^ 


S 2 (k) ; 

k^ < k < k^ 


S 3 (k) ; 

k^ < k < k^ 

(2.4) 

S 4 (k) ; 

k~ < k < k 
3 v 


S 5 (k) j 

k < k < oo 
V 



S (k) = « e -P(e 2 /(Ui9.5(u*)> V) 

1 2k J 

, 2 
k 2 u *m 

for 0 < k < k = (2.5) 


( 2 . 1 ) 
( 2 . 2 ) 

(2.3) 


S 2 (k) " “ 1/2. 5/2 

2k j k 

for k^ < k < k^ 


( 2 . 6 ) 
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for k^ < k < 


(2. 7) 


S 3 (k) 

S 4 (k) 

S 5 (k) 


aP < u ») 

2k 3 p k 3_p 


aD(u,; 

2k 3 

_ 3. 6 

Eu * k m 


vgk' 


for 


for 


k~ < k < k 

3 v 


k < k < oo 
v 


( 2 . 8 ) 

(2. 9) 


where S^(k^) = ^(k^) $ ^(k^) - S^k^) ; S^(k^) - S^(k^) and 

S.fk ) = S c (k ) ; and where a = 8.1 X 10 (3 = 0*74, E = 1.445 X 10 , 

4' v 5 V v' 

— =1.473 x 10 -4 , k, = k(6ir) = 0.359, k, = k(10u) S 0.942 

vg 2 ^ 

k = (gp/T)^ s 3.63 (t is the surface tension), 12 cm/sec, 

(2 

p r°g 10 ( k 3 /k 2 l 


k 

V 


2Eu , 

vg aD(u y 


1/6 


m 


.5756 

[D(u # )l 1/6 


( 2 . 11 ) 


Dfu^) = (1.274 + .0268 u^ 


-5 2 ^ 

+ 6.03 Xl0 uj 


( 2 . 12 ) 


and U in c (u ) is the wind at 19.5 meters for a friction velocity of 
19. 5 * 

in a neutral atmosphere. Ufu^) will often be used for U 19 5 ( u *)- 

These five ranges of definition can be substantiated by means of 
data obtained under a wide variety of conditions. S^(k) can be named 
the gravity wave-gravity equilibrium spectral range; S 2 (k) is the 
Kitaigorodskii spectral range; S 3 (k) can be called the Leyk in- Rosenberg 
spectral range; S 4 (k) is the capillary spectral range; and S 5 (k) is the 
Cox viscous cut-off region. 
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Anisotropy 


The definition of F(k, $) given by equation (2.3) was based on 
the work of Cote* et al. (I960) in project SWOP. It turns out to be inadequate 
for the higher wave number ranges, such as the capillary range. For 
these higher wave numbers not much was known then about the angular 
properties of the spectrum. 

For these high wave numbers, the exponential term in equation 
(2.3) is one, and the second term can be rewritten as equation (2.13). 

— [1 - ^ + A (cos $)^] = ^ [l + | cos 2$] = ~ [l + a^ (k) cos 2$] (2.13) 

In this form, it is seen that this term represents the first 
term in the even Fourier series representation of F(k, $) over 

TT TT 

- and that a more general representation would be given by 

equation (2.14). 


F(k f «) 


T T 


N 

1 + Fa (k; u ) cos 2n$ 
, n 
n= 1 


(2.14) 


Some results on the first two coefficients of this series 
a^(kj, u^) and a^(k, u^) can be obtained and will be given in a later 
section. 
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WIND PROFILE RELATIONSHIPS 


The average wind varies with elevation in the planetary boundary 
layer. Since different experimenters use anemometers at different ele- 
vations for free atmosphere measurements, and since data obtained in 
wind-water tunnels have to be interpreted for open ocean conditions, the 
procedures for interrelating the various measurements have to be given. 

The spectra are given in terms of u^ and g(u^) where u * 

is the friction velocity defined by = \fr / p , where t is the wind stress 
and where U, n r (u ,) is the wind measured at 19.5 meters in a neutrally 
stratified atmosphere. The equations relating these quantities according 
to Cardone (1969) are given by equations (3. 1) and (3. 2). 

U = rf ln(z/z ) (5. 1) 

z k o 

z = liM + 4.28 X 10" 5 u 2 - .0443 (3.2) 

o u , * 

* 

where z is the height of interest in centimeters (1950 in this case), 
k 1 is von Karman T s constant ( = 0.4) and U and u^ are expressed in 
centimeters / sec. 

Also since many measurements are referred to an elevation ot 
10 meters, the interrelationship between u^, U 10 ( u 5!< ) and Ujg 
is needed. Table 3.1 gives the values of these quantities, and Figure 3.1 
is a graph of g and U 1Q as a function of u^ for a neutral atmosphere. 

The range from 5 cm/sec to 138 cm/sec corresponds to wind speeds 
at 19.5 meters of 1.25 m/s (2.4 knots) to 27.0 m/s (52.4 knots). 

The wind at these two anemometer heights would be greater than 
the amount tabulated for a stable atmosphere to produce the same u s , c at 
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Table 3.1 Wind speed at 19.5 meters and 10 meters in 
meters/sec as a function of in cm/sec 
for a neutral atmosphere. 


% ' 

U 10 

U 1 9.5 

% 

U 10 

U 1 9.5 

5 

1.16 

1.24 

72 

15.45 

16.65 

6 

1.43 

1.53 

75 

15.92 

17.17 

10 

2.62 

2.78 

78 

16.38 

17.68 

12 

3.26 

3.46 

80 

16.69 

18.02 

15 

4.28 

4.53 

84 

17.29 

18.69 

18 

5.30 

5.61 

85 

17.44 

18.86 

20 

5.93 

6.27 

90 

18.18 

19.68 

24 

6.98 

7.38 

95 

18.91 

20.49 

25 

7.21 

7.63 

96 

19.05 

20.65 

30 

8.24 

8.74 

100 

19.62 

21.29 

35 

9.18 

9.77 

102 

19.90 

21.61 

36 

9.37 

9.97 

105 

20.32 

22.08 

40 

10.10 

10.75 

108 

20.74 

22.54 

42 

10.45 

11.15 

110 

21.01 

22.85 

45 

10.98 

11.73 

114 

21.56 

23.47 

48 

11.51 

12.31 

115 

21.70 

23.62 

50 

11.85 

12.69 

120 

22.37 

24.37 

54 

12.53 

13.43 

125 

23.03 

25.12 

55 

12.70 

13.62 

126 

23.15 

25.27 

60 

13.53 

14.53 

130 

23.68 

25.85 

65 

14.34 

15.43 

132 

23.54 

26.15 

66 

14.50 

15.62 

135 

24.33 

26.58 

70 

15.14 

16.31 

138 

24.71 

27.01 
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Fig. 3.1 U 10 (u + ) and Ujq gfu^) as a function of u^ in a neutrally- 

stratified atmosphere. 
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the surface and less than the amount tabulated for an unstable atmosphere. 
If possible, where needed, corrections for the effects of atmospheric 
stability have been made in interpreting the data. 


ZO 


THE GRAVITY WAVE-GRAVITY EQUILIBRIUM SPECTRAL RANGE 
Source 

The spectral range defined by equation (2.5) is obtained from 
the results of Moskowitz (1964), Pierson and Moskowitz (1964) and 
Pierson (1964) by transforming the frequency spectra given by them to 
wave number spectra. U = U^ ^(u^) represents the wind at 19.5 meters 
as described above. 

The equilibrium constant 

The equilibrium constant is usually found by attempting to fit 
the high frequency portion of a gravity wave spectrum to a curve with a 
slope of -5 on a double log plot under the assumption that 

S(co) = ag^/oo^ (4. 1) 

in this range. There has been some discussion concerning the possib- 
ility that the constant a in the Phillips equilibrium gravity wave range is 
fetch dependent as, for example, in Longuet-Higgins (1969), and Mitsuyasu 

(1969). For open ocean conditions, the preponderance of the data now 

-3 

suggests that a is a constant very nearly given by 8.1 X 10 as ob- 
tained by Pierson and Moskowitz (1964). 

Longuet-Higgins et al. (1963) gave a value of a equal to 

_ 3 

8x10 for open ocean conditions and a large fetch. 

3 

Ewing (1969) found that a = 8 X 10 fitted spectra obtained 
from a series of observations taken with the NIO buoy. 

Schule, Simpson, and DeLeonibus (1971) in a study of laser 
altimeter wave records in wave number space found that a equalled 
8.4 x 10" . The value of a does not seem to depend on fetch because the 
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fetch varied over one order of magnitude in their measurements. 


More recently, DeLeonibus et al. (1973) have summarized re- 
sults obtained at Argus Island and by using a bow-mounted wave height 

sensor in the Norwegian Sea. For Argus Island data, a was found to 

-3 - 3 

equal (7.8 ± 1.6) X 10 and for the Norwegian Sea a value of 8.3 x 10 

was found. 

A recent study by Saenger (1972) also supports a value of a 

-3 -3 

equal to 8.1 x 10 , or 8.0 x 10 , with one important exception to be 

described below. The wave recording instrument used was a heave ac- 
celeration, pitch and roll buoy with extremely carefully designed sensing 
elements. The high quality data obtained yield S(f), or equivalently, S(oo), 
from 0.5 hertz (f) to 1.0 hertz with a small roll-off in frequency response 
beyond about 0.7 hertz. Spectra were obtained for winds of 3, 6, 6.5, 

9.5, 11, 13,5, 14.5, 15, 17 and 18 knots. Over the frequency range from 

0.35 to 0.7 hertz, and for winds of 9.5 knots, or more, the spectra co- 

-3 

incide within sampling variability and correspond to a = 8.1 x 10 

within the accuracy of the plots. The spectra for the waves when the winds 

-3 

were 6 knots and 3 knots were well below this value with a = 3.2 x 10 

_ 3 

for 6 knots and 1.6 x 10 for 3 knots. 

This result is important because something even more dramatic 
happens to capillary waves at a wind speed of about 3.5 m/s (7 knots). The 
Phillips } equilibrium gravity wave range may only exist if the winds exceed 
3.5 m/s in a neutral atmosphere. 

In summary, for winds greater than about 3.5 m/s, Pierson 

_ 3 

and Moskowitz (1969) gave 8.1 x 10 for a, Longuet-Higgens et al. (1963) 
gave 8 x 10‘ 3 , Ewing (1969) gave 8 x 10 3 , Schule et al. (1971) gave 
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8,4 X 10” 3 (or Its equivalent in wave number space), DeLeombus et al. 

(1973) gave (7.8 ± 1.6) X 10 -3 and 8.3 X 10 3 , and data from Saenger 

-3 

(1972) are consistent with 8.1 X 10 ! 

A value of 8.1 X 10” 3 is used here with a variation from 
7.8 to 8.4 X 10 -3 , implying a 4 percent change in spectral values. If the 
spectral shapes were unchanged, this would imply a completely undetect- 
able (due to sampling variability) change in significant wave height of plus 
or minus 2 percent if the upper or lower ranges of this value were used. 

The SWOP project (Cote" et al, I960) is often cited for a value of 
a equal to 13.3 X 10 -3 . This value depends on the adequate calibration of 
a wave buoy that moved up and down in the water and may perhaps be in 
error for this reason. 

Other values of a often cited are 14*8 X 10 for a 70 meter 
fetch due to Burling (1959), 12.1 X 10” 3 for a 55 meter fetch due to 
Hicks (I960), and 10.4 X 10" 3 for a 2.5 km fetch due to Kinsman (I960). 
An explanation for these higher values will be given later m this study. 

Relationship to wave forecasting methods 

This portion of the spectrum is based on a selected subset of 
many hundreds of ocean wave records obtained by the Tucker shipborne 
wave recorder (Tucker, 1956) as analyzed by Moskowitz, Pierson and 
Mehr (1962, 1963). Fully developed seas with a spectrum given by this 
form are not frequent over the open ocean. This form, however, is the 
basis for the numerical spectral wave forecasting model described by 
Pierson, Tick and Baer (1966) and more recently by Pierson (1971) and 
Pierson and Cardone (1971). Papers by Cardone (1969), and Adamo, 
Baer and Ho s me r (1968) describe features of this hindcasting and fore- 
casting model. 
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An earlier version of this model has been tested on simulated fore- 
casts with results described by Bunting (1970) and Bunting and Moskowitz 
(1970). Both observations and these simulated forecasts show that the 
wave spectra depart considerably from the spectrum that would be present 
as calculated from the instantaneous local wind. 

Tests have shown, however, that this portion of the wave spectrum 
can be both forecasted and hindcasted with considerable accuracy on the 
basis of the wind field as it varies from hour to hour over the oceans. 

If the low wave number, low frequency part of the wind wave and swell 
spectrum had any effect on the amount of backs cattered radiation in radar 
measurements, then these forecasts and hindcasts can provide the required 
information on this part of the spectrum. 
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THE CAPILLARY SPECTRAL RANGE 


Data sources 

The capillary spectral range given by equation (2.8) has been 
determined from data obtained in three different wind water tunnels for 
a variety of fetches. A range of friction velocities from 3.3 to 
170 cm/sec corresponding to winds over the ocean varying from nearly 
calm to more than 29 m/s were covered by the data. 

The three data sources consisted of a set of wave records in 
digitized form obtained at New York University, a second set of records 
in graphical form provided by Dr. Yoshiaki Toba of the Geophysical Insti- 
tute of Tohoku University, which were originally obtained at Kyoto Uni- 
versity, and a third set of digitized records obtained by Dr. A. J. Suther- 
land and provided by the Department of Civil Engineering of Stanford 
University. The first set were for fetches of 2.3 to 6.4 meters and 
friction velocities of 3.3 to 19 cm/sec. The second set were for fetches 
of 4.5 to 13.6 meters and friction velocities of 33 to 115 cm/sec. The 
third set were for fetches of 1.1 to 18 meters and friction velocities of 
45 to 170 cm/sec, with data actually used having been for fetches of 
8.7 to 18 meters. The Nyquist frequencies for the three sets of data were 
78, 125 and 125 radians /sec respectively. 

A critical friction velocity 

The New York University data covered a range of friction velo- 
cities from nearly calm to 19 cm/sec. For the lower wind speeds, the 
wind hardly ruffled the surface, but slight disturbances could be detected 
that were capable of being recorded on a high gain oscillograph. At a 
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critical friction velocity of about 12 cm/sec the waves grew rather suddenly, 
and the data on their properties merge with the data from the other two 
sources as will be shown. 

The problem remains though as to what the spectrum would be of a 
fully developed sea for winds corresponding to friction velocities less 
than 12 cm/sec. Our results will be definitive only for u > 12 cm/sec. 
Perhaps winds less than this can only ruffle the surface and equation (5.1) 
given by Pierson and Moskowitz (1964), 

H = 0.0212 (5.1) 

where is in meters and Uj 9 5 is in meters/sec only holds for 

3 

winds greater than about 3.5 m/s at 19.5 meters, corresponding to a signi- 
ficant wave height of 0.25 meters or 11 inches. 

For wind speeds less than the wind that corresponds to a friction 
velocity of 12 cm/sec, the gravity wave spectrum may be very low under 
the idealization that a constant wind has blown for a long time over a large 
fetch. The minute undulations that were measured for these light winds 
yielded spectra for the capillary range that corresponded to disturbances 
one one -hundredth the height of the waves that were measured after the 
friction velocity exceeded 12 cm/sec. The data from Saenger (1972) could 
lead to the conclusion that the waves for a fully developed sea would be 
5 inches high instead of 8 inches for a 6 knot wind and only about 1 inch 
high in the gravity range for a 3 knot wind . There will, of course, 
be swell present. 
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A capillary equilibrium range? 

Phillips (1966) has discussed the pos sibility of two different equi- 
librium ranges : one for gravity waves and the other for capillary waves* 

In the previous section of this paper, a possible equilibrium range in 
the gravity wave portion of this spectrum was shown with a constant equal 
to 8.1 X 10 for winds above 3.5 m/s. 

Phillips (1966) wrote that the capillary spectrum should be of 
the form B*/iT 4 , or in polar coordinates after integrating out the angle 
variable, B/k^, where B is a constant. Phillips (1966), page 11 "7, can 
be quoted as follows: 

"It should be noted that an infinitesimal capillary wave 
cannot be in a state of stable equilibrium in which the energy 
flux from the wind balances the viscous dissipation. For ac- 
cording to the linear theory, these two energy fluxes are in 
fixed proportion, each being proportional to the square of the 
wave amplitude. A disturbance either dies away or grows until 
it is limited by one or the other of these finite amplitude effects." 

The data described in this section establish that there is no equi- 
librium range in the capillary spectrum as computed from wave records 
obtained in a wind-water tunnel. This fact was recognized by Sutherland 
(1967) who wrote that "the present data do not show a true equilibrium 
range as predicted by Phillips (1958)." 

The capillary spectrum is a function of wind speed. An adequate 
theory for this portion of the spectrum of a wind roughened sea is lacking. 
The suggestion given by Pierson, Jackson, Stacy and Mehr (1971) is not 
adequate because spectral components with the least viscous effect grow 
the most slowly. 
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Data analysis 

If the capillary spectrum is of the form given by equation (5. 2) 

S(k) dk = dk 
k 


(5. 2) 


and if there are no currents* the use of equation (1. 13) can transform the 
wave number spectrum to a frequency spectrum as in equation (5. 3). 


S(w) dw = 5__ 


(kM) 


3 dw 


dw 


(5.3) 


Let the constant B be assumed to be equal to 4.05 X 10~^,and let 
the spectra estimated from the records described above be designated 
by S(w; u^,F). If then the capillary spectrum is of the form given by 
equation (5.3) the function of w given by equation (5.4) ought to equal 
some constant over some portion of the capillary range of frequencies 
in the spectrum. 

.3 


H(w; u^, F) = 


S(w,u,, F)(k(w))' 


(5.4) 


(a / 2) (dk(w) / dw) 

The function, H(w; u s! _, F) was computed and graphed for a wide 
range of estimated spectra S(w; u ,, F), for the ranges of u and F 
given above. H(w; u^, F), in general, was small for low values of w, 
increased sharply to some constant value at some Wj , and leveled off 
to nearly a constant value after that w. To determine this constant value, 
the quantity D was computed as in equation (5. 5). 


w 


D(u.F) = 


1 


N 


w M - u, 

N i w 


/ H(w; u , F) dw 


(5.5) 
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The value of varied as a function of both and F . If F 2 >Fj 

and u^ was constant, then C0 ^2 <00 11* ^ an ^ ^ 

Fj = F 2 , then <^ 12 < l • The spectra of the waves In the wind water 
tunnel grew toward lower frequency with both Increasing fetch for the same 
friction velocity and with Increasing friction velocity for the same fetch. 

Selected graphs of H(oo; u^, F) are shown from the Sutherland data 
in Figure 5.1 for different fetches and values of u ^ . The values of 

and gj , which is where the effect of digitization noise resulted in the upper 

N 

cut-off value, are shown along with the resulting value of D(u >;c ). The plotted 
points as connected by a zigzag line are the result of 21 -point non-overlapping 
averages of the Fast Fourier transform of the data. The 90% confidence 
interval is shown, and the variability of H(cj; u^, F) is consistent with 
this confidence interval. The value of D(u^) is insensitive to the exact 
choice of . 

A total of 67 spectra were obtained for varying friction velocities and 
fetches, 27 from N. Y. U. , 23 from Dr. Toba, and 17 from Dr. Sutherland 
(1967). These spectra were processed to obtain H(a>; u^, F) and D(u^, F). 

2 1 i 

The values of u , F, ( u^ /gF) 2 , D, and D 2 are shown in Table 5.1. 

Sutherland (1967) stated that the spectra he obtained were not good 
at high frequencies. This point was overlooked in the initial analysis of 
the data so that the values of D as plotted by Pierson, Jackson, Stacy and 
Mehr (1971) came out, without a needed correction, to be higher than the 
values obtained from the Toba data. 

It was found that the Sutherland data had white noise at high fre- 
quencies, probably caused by digitization effects. This white noise was 
easily removed by finding its magnitude at the high frequencies and 
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Fig. 5.1 Selected graphs of H(wj u^, F) from Sutherland data. 
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Table 5.1 


2 — — 

Values of u^ (cm/sec), F (meters), (u^/gF) 2 , D, and D 2 
for the N. Y. U. , Toba, and Sutherland data sets. 



Data set 

u . 
❖ 

F 

(u */ gFF 

D 

D* 

1 . 

NYU 

3.3 

2.3 

6.9 X 10" 3 

3.8 X 10" 5 

-- 

2 . 

1 ! 

3.3 

4.8 

4.8 x 10' 3 

4.5 X 10" 5 

-- 

3 . 

11 

3.3 

6.4 

4.2 x 10“ 3 

7.2 x 10" 5 

-- 

4 . 

1! 

4.7 

2.3 

9.9 x 10' 3 

9.6 x 10" 5 


5 . 

1 1 

4.7 

4.8 

6.9 x 10" 3 

6.9 X 10" 5 

-- 

6 . 

1 1 

4.7 

6.4 

5.9 x 10' 3 

8.5 x 10" 5 

-- 

7 . 

II 

5.9 

2.3 

1.2 x 10 -2 

2.4 x 10" 4 

-- 

8 . 

1 1 

5.9 

4.8 

8.7 X 10 " 3 

9.5 x 10" 5 

-- 

9 . 

If 

5.9 

6.4 

7.5 x 10 " 3 

1.35 x 10" 4 


10 . 

H 

7.4 

4.8 

1.1 X 10' 2 

7.75 X 10" 4 

-- 

11 . 

1 1 

7.4 

4.8 

1.1 X 10" 2 

1.22 x 10" 4 

-- 

12 . 

i r 

7.4 

6.4 

9.4 X 10 " 3 

1.38 X 10" 4 

-- 

13 . 

i ? 

9.2 

4.8 

1.36 X 10" 2 

4.32 x 10" 4 

-- 

14 . 

1 1 

9.2 

6.4 

1.15 X 10" 2 

1.42 X 10" 4 

-- 

15 . 

1 1 

9.2 

6.4 

1.15 X 10" 2 

1.26 X 10" 4 

-- 

16 . 

i f 

9.9 

6.4 

1.25 X 10" 2 

3.14 x 10" 4 

-- 

17 . 

1 1 

9.8 

6.4 

1.25 X 10" 2 

2.68 X 10" 4 

-- 

18 . 

i i 

12.1 

4.8 

1.79 x 10" 2 

7.8 X 10" 3 


19 . 

1 ! 

12.1 

6.4 

1.53 X 10' 2 

1.56 X 10" 1 


20 . 

1 i 

12.3 

6,4 

1.55 X 10" 2 

2.78 

1.67 

21 . 

1 ' 

13.3 

4.8 

1.94 x 10" 2 

2.76 

1.66 

22 . 

1 T 

13.8 

6.4 

1,74 X 10" 2 

3.08 

1.75 


[continued] 
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Table 5.1 (cont. ) 



Data set 

u * 

F 

K 2 /gF ) 2 

D 

l 

D3 

23 . 

NYU 

14.2 

4.8 

2.07 X 10" 2 

2.50 

1.58 

24 . 

M 

14.2 

6.4 

1.79 X 10' 2 

3.20 

1.79 

25 . 

i i 

16.2 

2.3 

3.41 x 10 -2 

1.69 

1.30 

26 . 

1 1 

16.2 

4.8 

2.36 X 10“ 2 

3.48 

1.86 

27 . 

ii 

16.2 

6.4 

2.04 x 10" 2 

3.64 

1.91 

28 . 

n 

18.9 

2.3 

3.98 X 10“ 2 

1.19 

1.09 

29 . 

1 1 

18.9 

6.4 

2.39 X 10" 2 

5.68 

2.42 

30 . 

Toba 

33.4 

4.5 

5.03 x 10 -2 

4.40 

2.10 

31 . 

1 1 

33.5 

10.0 

3.38 X 10 -2 

4.19 

2.05 

32 . 

r i 

36.7 

13.6 

3.18 x 10‘ 2 

4.48 

2.12 

33 . 

1 1 

38.9 

13.6 

3.37 X 10" 2 

4.00 

2.00 

34 . 

1 1 

45.2 

6.9 

5.49 X 10" 2 

6.50 

2.55 

35 . 

1 1 

48.9 

13.6 

4.23 x 10" 2 

7.58 

2.75 

36 . 

f i 

49.3 

6.9 

5.99 x 10" 2 

9.76 

3.12 

37 . 

1 1 

56.5 

6 . 9 

6.87 X 10“ 2 

7.15 

2.67 

38 . 

1 1 

61.5 

13 . 6 

5.33 X 10" 2 

7.15 

2.67 

39 . 

1 1 

63.2 

6 . 9 

7.68 X 10" 2 

9.68 

3.1 1 

40 . 

1 1 

64.8 

10 . 0 

6.54 x 10" 2 

12.16 

3.49 

41 . 

1 1 

73.7 

13 . 6 

6.38 x 10“ 2 

13.68 

3.70 

42 . 

r ! 

77.4 

4 . 5 

1.17 X 10" 1 

15.5 

3.94 

43 . 

1 1 

78.6 

6 . 9 

9.55 X 10' 2 

16.56 

4.07 

44. 

1 1 

88.1 

6 . 9 

1.07 x 10' 1 

27.41 

5.24 


[continued] 
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Table 5 

1 (cont 

) 



Data set 

u * 

F 

(u^/gF)* 

D 

1 

D 2 

45. 

Toba 

92.6 

13.6 

8.02 x 

io" 2 

14.9 

3.68 

46. 

H 

102.0 

4.5 

1.54 X 

10- 1 

29.10 

5.39 

47. 

1 1 

104.0 

10.0 

1.05 X 

10- 1 

23.90 

4.89 

48. 

it 

112.0 

10.0 

1.13 x 

io- 1 

32.4 

5.69 

49. 

H 

114.0 

13.6 

9.87 X 

io- 1 

21.7 

4.66 

50. 

1 1 

115.0 

4.5 

1.73 x 

io- 1 

24.2 

4.92 

51. 

Suthe rland 

37.2 

18.0 

2.8 x 

io- 2 

7.44 

2.73 

52. 

u 

40.3 

11.9 

3.73 x 

io- 2 

6.20 

2.49 

53. 

1 1 

40.1 

10.3 

3.99 x 

io- 2 

9.25 

3.04 

54. 

r i 

40.3 

14.8 

3.35 X 

io- 2 

6.57 

2.56 

55. 

1 1 

52.4 

17.9 

3.95 x 

io- 2 

6.30 

2.51 

56. 

1 1 

56.5 

10.3 

5.62 x 

io- 2 

7.98 

2.82 

57. 

! 1 

56.7 

11.8 

5.27 X 

io- 2 

6.70 

2.59 

58. 

1 f 

56.7 

14.9 

4.69 X 

io -2 

6.60 

2.57 

59. 

1 1 

84.2 

10.3 

8,38 x 

io " 2 

11.8 

3.44 

60. 

! 1 

84.8 

11.8 

7.88 X 

io - 2 

12.60 

3.55 

61. 

n 

85.8 

14.9 

7.10 x 

io - 2 

12.58 

3.55 

62. 

r i 

86.2 

17.9 

6.51 x 

io - 2 

11.43 

3.38 

63, 

r t 

102.4 

8.7 

1.11 X 

io - 1 

22.07 

4.70 

64. 

1 1 

104.7 

10.3 

1.04 x 

10" 1 

24.21 

4.92 

65. 

I ! 

105.8 

1 1.8 

9.84 x 

io - 2 

23.46 

4.84 

66. 

ft 

169.4 

8.6 

1.84 x 

io " 2 

49.47 

7.06 

67. 

1 1 

170.2 

10.3 

1.69 X 

io" 1 

60.96 

7.81 
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subtracting it from the spectrum before calculating H(oj; u^, F) and D(u^). 

The Toba data were provided to us in the form of strip charts 
of wave height as a function of time. These charts were magnified by 
projecting them on a wall across the room and the data were read off 
very accurately. There was no need for a correction for white noise 
digitization errors in the Toba data processed in this way. 

The variation of D(u,) as a function of u is shown on a double 
log plot in Figure 5.2. There are three anomalous points for D greater 

than one in the NYU data. For u, less than 12 cm/sec, D increases 

❖ 

-4 -3 

from values near 5x10 to values near 4x10 at u^ equal to 
10 cm/sec. Near 12 cm/sec, values of D scatter over four orders of 
magnitude, and for u^ just greater than 12 cm/sec the values of D 

*i c 

cluster around three. Since the wind at anemometer level is not a very 

good measure of u^ when the winds are light, because atmospheric 

stability plays an important role* this plot suggests a reason for those 

radar measurements that found strong effects of varying wind speed on 

radar sea return for light winds. The coded points attributed to Bradley 

(1971) will be discussed in a later section. 

There is a data gap between 20 and 30 cm/sec. Above 30 

cm/sec the Toba and Sutherland data would be indistinguishable were 

it not for the coded points. The data show a constant increase in D 

from 4 to 60 as the friction velocity increases from 33 to 170 cm/sec. 

Although the frequency at which H(oji u >t< , F) first levels off to a 

constant value is fetch dependent, the value of D does not appear to be 

fetch dependent. This is shown by Figure 5.3 in which D is plotted 
2 — 

against (u^ /gF) 2 on roughly the same horizontal spacing. The scatter 
in the points is increased compared to Figure 5.2, whereas if fetch were 
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truly involved, the scatter should have decreased. The square root of the 
reciprocal of the dimensionless fetch varies over nearly two orders of 
magnitude. One could interpret this plot to infer that D is proportional 
to ((u^/gF)" 2 )* 1 for ( u JV g F ) 2 greater than 1.2 X l(f 2 . Only the in- 
creased scatter when compared with figure 5.2 leads to the preference to a 
wind speed dependence only. Some high quality measurements of waves 

for a fetch of one kilometer or so for the same frequency range could 

* 

settle the question completely. Note that the coded points attributed to 
Bradley (1971) in Figure 5.2 are not shown on this figure. 

A relationship between D and u is needed to provide an analyti- 
cal representation for the spectrum. It should be recalled that the pri- 
mary measurement was wave elevation as a function of time, and that 
spectra are related to squares of wave elevation. Since the spectrum for 
a long enough fetch is more or less defined over a constant wave number 
range the square root of D would be a measure of the wave height in 
this band, 

1 

Figure 5.4 shows a plot of (D(u^)) 2 versus u^ on linear scales. 

The significant height of the capillary waves can be given approximately 

by equation (5.6) because the wave number range for the capillary portion 

of the spectrum extends from 0.92 radians /cm to about 10 radians /cm 

on the basis of results to be given later. An upper limit of the integral 

ranging from 6 to 12 has only a small effect. 

r 10 4.05 x 10" 3 D(U ) dk A 

H m 34 / = 0. 1 95 (D(u )) 2 (5.6) 

1/3 -0.94 k 3 _ * 

1 

A second vertical scale shows that the variation of (D(u^)) 2 from 1 to 8 

>;< 

See page 81. 
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corresponds to a variation in significant wave height^ from 0,3 to 1,6 cm. 

The least squares regression of [D(u^)] 2 on is shown in 
Figure 5.4 for > 12 cm/sec. 

The defining equation for [Dfu^)] 2 is~ given on the figure and 
by equation (5.7) (or equation (2,12)) for u > 12 cm/oec. The three out- 
rider points in the NYU data are included. Their omission would not change 
the equation very much. 

(D(u + ))1 = 1.247 + 0.0268(u j{c ) + 6.03 x l(f 5 (u*) 2 (5.7) 

This form of presentation also shows something about the size of the 
scatter. The significant capillary wave height scatters by about plus or 
minus 0.1 cm about the regression line over the full range of u^ from 
12 to 170 cm/sec. 

There would have been other ways to obtain an expression for 
D(u^) as a function of u^, but the way chosen appears to fit D(u^) as 
given in equation (5.7) versus u^ on both linear and double log scales al- 
most as well as if the regression curves had been derived to fit the data 
on these scales as shown in Figures 5.5 and 5.6. Note the coded points 

attributed to Bradley (1971) in Figure 5.6, which will be discussed later. 

1 

Given the relationship between u^ and ^ , D 2 can be plotted 

i 

against U J9 ^ as in Figure 5.7. D 2 varies from about 1.5 at 3.5 m/s 
to 7 at 30 m/s and the wave height varies from 0.3 to 1,4 cm. A wave 
6.8 cm long and 1,4 cm high is shown as an insert. Figure 5.7 also 
shows the sharp increase in wave height near 3.5 m/s. Slightly below 

^The significant wave height is defined to be the average of the heights 
of one-third highest waves in a wave record, in this case after the other 
wave number ranges have been filtered out. 
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this speed, the capillary wave height decreases sharply from 0.3 cm 
to 0.003 cm for this range of wave numbers. 

Ph otographs of capillary waves 

Figures 5.8, 5.9 and 5.10 show 32 cm long photographs of the 
waves generated at Kyoto as furnished by Dr. Xoba for nominal wind tunnel 
speeds varying from 5.1 to 12.1 m/s and for fetches of 4.5, 6.9, and 
10.0 meters. The high wave number, short wave, structure of the surface 
becomes rougher with increasing wind speed at all three fetches, and 
this short wave structure is nearly the same for a given wind speed at 
all three fetches. The figures also show the crossing capillary waves in the 
proper locations as described theoretically by Longuet-Higgins (1963) 
for the lower wind speeds, and as predicted by Phillips (1957) in the reson- 
ant theory of wave generation. These crossing waves at low winds will 
play an interesting role in providing a refined version of equation (2.3) 
in a later section. The surface becomes uniformly rough for the 
higher speeds. It is difficult to imagine any reason why these scales 
of capillary waves should become lower for open ocean conditions. 
Moreover, these photographs support the procedures used in going back 
and forth from frequency to wave number space. The wave number 
spectrum clearly grows with increasing wind speed for waves between 
6.8 and 0.63 cm in length. 

Catspaws 

One of the interesting features of these results, in particular 
Figure 5.2, is that they can be used to explain ’’catspaws M . Catspaws 
are patches of roughened water scattered over a relatively smooth surface. 
They form, move along for a short distance, and disappear to be replaced 
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by other "catspaws" somewhere else over the water. If the wind over the 
water is describable as ’’light air”, the friction velocity at the surface 
will be less than 12 cm/sec over nearly the whole surface. Turbulent 
gusts will produce small isolated patches over the water where u^ 
will exceed (12 cm/sec). The spectrum of the capillary waves in 

these patches will be four orders of magnitude higher than the spectrum 
of the low undulations in the nearby areas. The waves in the "catspaws” 
will be 100 times higher than the waves outside the catspaws. 

Wave number spectra in the capillary range 

Dr, Alan Schooley provided us with photographs for two different 
wind speeds of the wavy surface along the glass side of a small wind- 
water tank. The fetch was very short. The elevation as a function of dis- 
tance along the glass side was read off and two wave number spectra were 

_3 

computed. The spectra varied as (£*) (see equation (1.16) and the dis- 
cussion in the section on the Leykin-Rosenberg range) with the one for the 
higher wind speed being considerably higher than the one for the lower 
wind speed at all values of £*, Since u^ was not known, exact quanti- 
tative agreement with the preceding equations could not be obtained, but 
the spectra behaved as the theory would predict and the photographs in 
this section show they should. 
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THE COX VISCOUS CUTOFF 


Cox (1958) measured the upwind -downwind slope of capillary 
waves as a function of time at a fixed point in a small wind-water tunnel 
with a fetch of 2. 14 meters* He was able to determine the spectrum of 
the slope of the waves out to very high frequencies. For the lowest 
wind speed for which data were obtained there was a sharp falloff of the 
slope spectrum by three orders of magnitude from ZOO radians /sec to 
800 radians /sec. For the next highest wind speed the three orders of 
magnitude falloff occurred between frequencies of 800 radians /sec and 
2000 radians /sec. 

From the available information, which consisted of a large size 
copy of the original plot of the data sent to us by Dr. Cox, it was possible 
to fit a function to the sharp cutoff region of these data. It was found that 
equation (6. 1) described this portion of the spectrum. The constant A 
equals 0.79Z and is not dimensionless. 


A u 


fS s {£) =_ x 


( 6 . 1 ) 


The frequencies, go > ZOO, were so high that the capillary 
relationships for wave number and frequency could be used. Equation (6.1) 
can be put in the form given by equation (6.Z) where A* is dimensionless. 


S (f) 
s ' ' 


A i 2 3 

A v k go u 
mm* 

7 

v f 


( 6 . 2 ) 


3/2 i 

Then by means of go = Zuf, go = — 5 k , and go = (Zgk ) z , 
7 2 m & m 


m 


it is possible to find k S(k) and finally S(k). The result is equation (6.3), 
(or equation (2.9)). 
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(6.3) 


u 3 k 6 

_ ,, . , n ,4,* * m 

S 5 (k) = 3(2 tt) A — 

k 



3 

The dimensionless form /vg has been constructed with v 
as one of the parameters because it yields a dimensionless result; ob- 
viously other combinations of k , g, v, and tu could also have been used. 

It would be more desirable to obtain both 9r|(t)/8x and 8r](t)/8y 

and proceed from it to S (k), but such data are not presently available so 

5 

that the procedure used will have to suffice. As far as the viscous cut- 

Z 2 

off region is concerned, the function k (cos$) S(k)F(k, $) has effectively been 

2 

treated as if it were k S(k)F(k, $). 

If S^(k^) is set equal to S 5 (k^), the wave number at which the 
k behavior shifts to the k behavior can be found. The value of k^ is 
given by equation (6.4). 


k 

v 


2Eu 3 

* 


1/6 


k 

m 


(6.4) 


Over the range of u^ corresponding to realistic winds increases 

slowly with increasing u^ . The value of k^ is typically near 10 and 
actually varies from 6 to 12, or so, as u varies from 12 to 48. 

The wavelength that corresponds to k = 10 is 0^6 3 cm which yields the 
upper limit of the integral in equation (5.6). 

Cox (1958) also found the equivalent of the threshold u at 
about 12 cm/sec described above. The change in boundary layer 
characteristics at the input end of the wind water tunnel seems to have 
produced a slight fetch effect. The 10 meter wind at which the waves 
were observed to increase suddenly, with appropriate remarks about 
the unreliability of the calculation, was found to be given by 2.9, 2.5, 
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and 2. 3 m/s as the fetch varied from 1. 47 through 2. 14 to 3. 12 meters. 
These values would correspond to values of u between 9 and 11 cm/sec 
if the estimates were correct. If they were more nearly 3. 5 m/sec, then 
the value of u^ would agree with the value deduced in the preceding analy- 
ses of the capillary range. 

The measurements made by Cox (1958) of the upwind -downwind 
slope as a function of time at a point were made for values of u A of 
284 cm/sec, 217 cm/sec, 144 cm/sec, and 31 cm/sec. These correspond 
to 19. 5 meter average winds of 45. 6 m/s, 37. 6 m/s, 27. 9m/s, and 
8. 98 m/s; and, in turn, the first two of these speeds corresponds to a 
wind in excess of hurricane speed and to a wind of nearly hurricane speed. 
The Sutherland (1967) data on elevation spectra for longer fetches only 
extends to about 180 cm/sec for comparison. 

Figures 6. 1 and 6. 2 show graphs of fS g (f) versus f as found by 
Cox (1958) plotted along with the theoretical graphs for the same quantity 
for an infinite fetch for these four different values of u^. The agreement 
is poor for u^ = 284 cm/sec both because there is no particular reason 
why the equation for D(u^) should be valid for such an extreme wind and 
because the fetch is so short that the waves might not have been generated 
anyway. 

For u^ = 217 cm/sec, the viscous cut-off range agrees fairly 
well but the observed spectrum is low compared to the computed one. 

For u yc = 144, the viscous cut-off range is correctly located and the 
spectra agree over portions of their range of definition within a factor 
of two. 
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Fig. 6.1 Theoretical and observed graphs of fS g (f) versus f for high 
values of u - 
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Fig, 6,2 Theoretical and observed graphs of fS (f) versus f for moderate 

values of u^ . S 

* 
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For u =31, the agreement Is really quite good over one 

>{< 

order of magnitude in f. The fetch is simply too short to generate those 
components with frequencies less than 6 hertz. 

The further investigation of the slope spectrum over a more 
modest range of wind speeds for longer fetches should provide useful 
refinements in the definition of the viscous cut-off range. Research 
on this problem by others is under way at the present time and the results 
can be used to improve on the analysis given here. 
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THE KITAIGORODSKII AND LE YKIN-ROSENBERG WAVE 

NUMBER RANGES 


The need to connect the gravity wa ve and capillary ranges 

From the preceding sections, a gravity wave spectrum has been 
obtained that has the property that 


S, (k) 

2k 


(7. 1) 


for some high wave number part of the gravity wave spectrum, 0<k<k^. 


A capillary spectrum given by 

a D(U * } 

S 4 (k) = ^ ~y~ 

k 


(7. 2) 


for k~ < k < k has also been found. 

■D V 

If u >; c * ^0, f° r example, corresponding to a 19*5 meter wind of 

12.7 m/s, then D(u^) equals 6.5 and S^(k) is 6. 5 times higher than 

the continuation of S^k) into the interval, k^ < k < k^ . Similarly, if 

u * = 100 cm/sec, corresponding to a 19.5 meter wind of 21 m/s, D(u ) 

equals 20 and S^(k) is 20 times higher than the continuation of Sj(k) 

into the interval, k 0 < k < k . 

3 v 

It thus becomes necessary to connect the Phillips equilibrium 
gravity wave range with this highly variable capillary range and to deter- 
mine the values for k^, k 2 and k^. Equations (2.6) and (2,7) represent 
the Kitaigorodskii range and the Leykin-Rosenberg connecting range 
respectively. 

The Kitaigorodskii range 

The range defined by equation (2.6) was first proposed by 
Kitaigorodskii (1961) who discussed an interval of small-scale isotropic 
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turbulence postulated to follow fin frequency) the Phillips equilibrium range. 

A quotation from this source is: 

"We advance the additional hypothesis that at a frequency 
somewhat above the equilibrium interval the small-scale 
turbulent motions are such that they are affected only negligibly 
by the force of gravity, g. A region of the corresponding scale 
of magnitude is still not affected by surface tension and 
molecular viscosity. Then 


S(w) = £(E max ,w) 


(7.3) 


whence 


-4 , 

S(w) ~ Vax“ W T~ 

CO 

"Equation (7.4) can serve as the essential particulari- 
zation of [the spectrum] at high frequencies (somewhat above 
the equilibrium interval). ” 

This .interval of isotropic turbulence grows like u^j whereas 
the capillary range grows like D(uJ as shown previously. Moreover, it is 
not yet clear where one interval begins and the other ends. 


The Leykin- Rosenberg results 

Leykin and Rosenberg (1970) seem to have answered the question of 
how to connect the gravity wave range to the capillary range by means of 
observations made in the Black Sea at a point 50 meters from the shore in 
a water depth of 15 meters for situations where fetch effects did not influence 
the shape of the spectra. Thirty-two wave spectra were obtained "of which 
19 were in absolute values and 13 in relative values". 

The data show a clear increase of the spectrum at 3, 5, and 7 hertz 
as the wind speed at a 20 meter anemometer height varied from 0 to 
11 meters/second. A plot of S(f) X f 4 with the area under S(f) normalized 
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to one from 1 to 7 hertz is shown for data grouped in wind speed ranges 
of 0 to 0.5 m/s, 0.5 to 2 m/s, 2 to 4 m/s, 4 to 8 m/s, and 8 to 11 m/s. 
Neither an f nor an f ' law is strictly correct over this full frequency 
range for the controlling parameters are wave number dependent and not 
frequency dependent and the frequencies span those for gravity waves and 

capillary waves. The low frequency end for low winds appears to follow 

-5 

an f law because of the absence of the Kitaigorodski range. All plots 
of the normalized function S(f) X f show a relative minimum at f= 3 hertz 
and all but one show a rn.axim.um. at f = 5 hertz with S(f) f^ increasing 
between 3 and 5 hertz. 

Beyond the 5 hertz peak the different spectra all still increase 
with wind speed roughly in the way the capillary range described above 
increased in the spectra studied from the wind-water tunnel. 

Leykin and Rosenberg (1970) make several pertinent observations 
about their results. One is quoted as follows: 

" . • » At wind velocities V > 0.5 m/sec, a spectral density 
peak appears in the spectra S j(f) at a frequency of 4 to 
6 Hz. With an increase in the wind velocity from 0,5 to 
11 m/sec the peak increases more than twofold, while its 
location on the frequency axis remains practically un- 
changed. . . . rf 

A second is quoted as follows: 

"In order to ascertain the uniqueness of the relation between 
S(f) and the wind velocity, we carried out several series of 
consecutive swells [wave] measurements; each series of measure- 
ments lasted several hours; during this time the parameters of 
the low-frequency swell (wave] components remained [invariant], 
while the wind velocity varied. The results of the measurements 
confirmed the decisive effect of the wind velocity on the variation 
of S(f). " 
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A third is quoted as follows: 

"The observed increase in the spectral density of S(f) with 
increasing wind velocity may be attributed either to amplifi- 
cation of nonlinear interactions in the low-frequency range 
of the spectrum or to direct transfer of wind energy to the high- 
frequency range [according to Phillips]. Since the frequency 
range considered was, as a rule, located far from the energy- 
carrying maximum (f Q = 0.1 to 0.2 Hz), and taking into account 
the decisive effect of the wind velocity on the variation of S(f), 
we may apparently give preference to the second of the above- 
mentioned mechanisms. " 

A final quotation is: 

"...All the results considered correspond to wind directions 
at which sufficient fetch paths were ensured for all the wave 
components; therefore the effect of the fetch path on the change 
in the shape of the spectra may be neglected. " 

These observations seem therefore to confirm the belief that 
the high wave number, high frequency portion of a wind roughened sea 
responds quickly, and rather uniquely, to variations of u^. The upper range 
of wind speeds covered by the data was only 8 to 1 1 meters /sec so that open 
water verification of the wind-water tunnel results for higher wind speeds 
is still required. There is a discrepancy between this paper and the wind 
tunnel data on the critical wind speed because changes are reported in the 
ranges of 0 to 0.5 m/s and 0.5 to 2 m/s that would not compare with 
wind tunnel results. 

The variations in the data, the minimum at 3 hertz, and the 
maximum at 5 hertz all suggest that the highest wave number for the inter- 
val of small-scale isotropic turbulence should be given by k = k(u>) =k (6 tt), 
and that the lowest wave number in the capillary range should be given by 
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k - k(w) - k (IOtt). These values correspond to wave numbers of about 
0.35 and 0.92 or wavelengths of about 18 cm and 6.8 cm. The interval 


between 0.36 and 0.94 can be called the Leykin-Rosenberg range. The value. 


0»94, for kj determines the lower limit of the integral in equation (5.6). 

For u* = 12, the capillary range is still higher than the Phillips 
equilibrium range in the gravity wave spectrum. The u^ dependent 
range should perhaps disappear so that kj = k £ . The Leykin-Rosenberg 
range could be made to connect the gravity range to the capillary range 
by a straight line on a double log plot of S(k) versus k. 

At k = k 2 , S 2 (k) should grow linearly with u^. From (7.4), 

Gg u 

S 2 (w) = T (7.5) 

CO 


and by means of oj 


2 


= gk one obtains 


s 2 (k) = 


Gu 

2g^ k 


* 

¥U 2 - 5 


(7.6) 


At k. = k_, and u* = 12 = u 


*m ' 

s 2 (k) = 


Gu 

*m _ or 

2?kf' 5 2k* 


(7.7) 


so that 


and S^k) has the form 


G 

u *m k ! 


U , 


(7.8) 


s 2 (k) = Tfj ^ 713 

2k ^2 *m k 


(7.9) 


This carve will intersect the curve defining S-^(k) at k - k^ , and 
imposing this condition, one finds that 
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Zkj J 


exp - 


P V(UK))^' 1 l%m) k i 2 ' 


(7.10) 


The exponential term can be neglected as typical values of the ex- 
ponent work out to be 10”^ or so, so that the term is essentially unity. It 
follows that 




(7.11) 


and that kj decreases with increasing wind speed as also given following 
equation (2.5), 

At k = k 2 = k (6ir), the value of S 2 (k) is now given by 


c ,, , a U * 1 

S Z k 2 = 2 — “T 
'to k^ 

and at k = k 3 = k (10ir), the value of S 4 (k) is given by equation (7.13). 

n D < u *> 

S 4< k 3> = 7 — T* 
k 3 


(7.12) 


(7.13) 


The Leykin-Rosenberg range will be defined by connecting the 
two points S 2 (k 2 ), k 2 , and S 4 (k 3 ), k 3 by a straight line on a double log 

plot. The result is 

lo %0 S 3 (k) " 1 °glO S 4 (k 3 ) _ 1 °glO S 2 (k 2 ) : 1 °glQ S 4 (k 3 ) 


log 1Q k - log 10 k 3 


lo ^10 k 2 “ lo ®10 k 3 


(7.14) 


This result can be written in the number of forms such as equation (7.15). 


.S,Wv l0S -° ^ 

log io(s 4 (it 3 )y ' k 2 

log 10 kT 


log. 


s 2 ( k 2 ) 

10 


)] 


(7.15) 
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and equation (7.16) 


S 3 (k) 

S^i^) 


= (k/k 3 ) P 


(7.16) 


where P = 


lo gl0(S z (k 2 )/ S 4(k 3 )) 

Tog^(k^7kp 


From the definitions of 


anc * ^ 4 ^ 3 ) ' e< l ua ll on (7.16) can be put into a final form given by 


S 3 (k) = § 


D(%) 


(7.17) 


where 

lo gl0( D < u *> / ( u * /u *m )) 

p - log 10 (k 3 /k 2 ) 


(7.18) 


f°r < k < . When k - k^ # equation (7.17) equals S^(k) and from equation 

(7.14), when k = k^, it equals S^k^). 


Comparison of Kitaigorodskii range with data from Kinsman 

Kinsman (I960) carried out extensive measurements of waves in 
Chesapeake Bay and documented the results thoroughly so as to show certain 
non -Gaussian features of the probability density function of wave elevation. 
The observations were of such high quality that Longuet-Higgins (1963) used 
the data to verify a theory he derived on the non-Gauss ian distribution 
of wave elevation. Kinsman (1961) reported on the results of his measure- 
ments and concluded that he could find no consistent evidence of an equi- 

-5 

librium range like oj . In fact, the average of all the spectra, normalized 

-4.5 

to f/f , came out like c o 
max 

The literature of that period for studies carried out in reservoirs, 
and other smaller bodies of water, with wave recording systems that could 
acquire data out to 1.5 Hz, or so, gave values such as from 
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Hicks (1960) and w“ 4 * 5 for Burling (1959). An equilibrium law according 
to agV was not clearly documented in these studies. Part of the dif- 
ficulty was that some of the spectra were normalized to f/f ax and 
that the Tujkey methods of spectral analysis smooth over a fairly large 
frequency range with a triangular weighting function. 

Four of the original data sets of wave elevation time histories 
as tabulated by Kinsman (I960) were augmented by adding on zeros at 
the end, computing a properly corrected Fast Fourier Transform, and 
finding non-overlapping averages so that the spectral estimates were indepen- 
dent for a square shouldered spectral window. The data sets were for 
records 017, 018, 093 and 094, with the first two corresponding to a 
u of 15.9 cm/sec and the second two to a u. of 41.2 cm/sec. These 
records were chosen because they corresponded to the lowest observed 
u^ and the highest observed u^ after correction for atmospheric stability 
effects according to the methods described by Cardone (1969). The pairs 
of spectra were averaged to obtain a single spectrum for the two different 

u , values. 

❖ 

These two spectra are shown on Figure 7.1 as labeled. The 
2-5 - 3 

heavy line shows ag go for a = 8. 1 X 10 . The dashed line shows S^(w) 

as given by equation (7.19), which is the transformation of equation (7.9) 
to frequency space, for u^ - 11.2. The capillary effects do not cause 
S ^(cjo) to depart from the theoretical gravity form. 

S 2 (w) = ag^u^/k^ (12)co 4 (7.19) 

The dashed line Is the Kitaigorodskii range for u^ = 41.2. The Kitaigorod- 
skii range for u^ = 15.9 cm/sec hardly shows up at all the lower right 
corner of the figure and is not graphed. 
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_ 5 

The equilibrium gj slope is well defined for u^ = 15.9 for the 
high frequency range, but the spectrum is low. Kinsman (I960) notes calibra- 
tion problems for Records 017 and 018 that could account for this. The 
= 41. 2 spectrum tracks the combined equilibrium and KLtaigo rodskii 
ranges within sampling variability, except at high frequencies where sampl- 
ing variability and aliasing clearly account for the difference between the 
computed spectrum and the theory. Normalizing the two spectra at f / f max 
and averaging would have masked the simple wind speed dependence postu- 
lated by Kitaigo rodskii and made it practically impossible to find. 

The data of Saenger (1972) used to document the value of a also 
support the Kitaigo rodskii range by showing its absence for frequencies less 
than 0. 7 hertz. For the wind speeds for which the observations were made, 
the Kitaigorodskii range should not have been observed. There was no evi- 
dence of wind speed dependence in the equilibrium range for wind speeds 
corresponding to u^ greater than 12. 

These results suggest that there is indeed an equilibrium range, 
but that both frequency boundaries shift toward lower frequencies with increas- 
ing wind. The Kitaigorodskii range may have been present in those observa- 
tions that reached high enough frequencies, but the data analysis procedures 
were not designed to find it. 

The ’’variable 1 ' constant for the equilibrium range thus appears to 
be the result of how much of a u >t< dependent contribution at high frequencies 
was averaged into the computation after perhaps an ^/^ max frequency shift, 
and the variable exponent appears to be dependent on how many spectra had 
an uT dependence and how many had an w dependence as masked by noise 
and sampling variability. The "constant" also appears to be fetch dependent 
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because shorter fetches imply smaller and more response wave record- 
ing systems that obtain data for higher wave frequencies to yield spectral 
values that include the Kitaigorodskii range, * 


Comparison of theory with the Leykin-Rosenberg data 

Leykin and Rosenberg (1970) averaged the spectra that were esti- 
mated by combining them into five groups according to different ranges of 
wind velocity. In turn, these five averaged spectra were normalized by 
dividing them by the variance to obtain 


S rel (f) = S(f) /cr 


(7.20) 


4 4 

This result was in turn multiplied by f , and S r ^(f) X f was plotted on 
a logarithmic scale versus f on a linear frequency scale and shown as 
Figure 5 in the cited paper. 

The values of S re ^(f) are quite sensitive to what happens at low 
frequencies and depend on such things as whether or not the point plotted at 

I hertz has a band width from 0.5 to 1.5 hertz or some other width. It is 
thus not possible, without the original data, to reproduce their data re- 
duction procedure by comparable calculations on the theoretical spectra. 

4 

The function, f S(f), was computed by means of appropriate trans- 
formations on the defining equations for a single wind speed within the range 
of the grouped winds of Leykin and Rosenberg (1970). The result was then 
multiplied by that factor that would make it agree with the corresponding 
graph of Leykin and Rosenberg at 3 hertz. 

The results for wind speed ranges of 2 to 4, 6 to 8, and 8 to 

II m/s as calculated for 3.5, 5.5 and 8.5 m/s are shown in Figure 7.2. 

The agreement for the 2 to 4 and 6 to 8 m/s is quite good with the theoreti- 
cal curves coming close to the observations over the range from 2 to 6 hertz. 


* R.A. Stacy will have additional results on this range in the near future. 
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Division of both curves by f would result in the value at f = 7 being about 
-4 

4 X 10 lower than the value at f = 1 so that even the curves for winds 
from 8 to 4 m/s are in fair agreement. 

Additional supporting data 

Quantitative data are difficult to obtain to support these two ranges. 
Much of the data has been processed in one way or another so as to make 
it difficult to go back to the basic observations and test these new spectral 
forms. Two sources can be cited, however, that give quantitative support 
for these ranges. One source is the figure given by Hess, Hidy and Plate 
(1969), and the other source is the graphs of the original spectra used to 
prepare the report by Dobson (1970). 

The figure in Hess, Hidy and Plate includes w = 6tt and go = 10tt 
in its range from co = 0.2 to go = 100. The scatter of points between go = 2tt 
and o> = 6tt shows a tendency to behave more like go ^ than go ^ as go = 6-rr 
ls approached from the left. The scatter of points between go = 6tt and 
go = 10t r illustrates the combined effect of wind speed variability in the 
Leykin-Rosenbe rg range and short fetches. Above go - 10tt, if the spectra 
obtained from Toba and Sutherland were used, they would show the strong 
wind speed dependence obtained in the section on the capillary range. 

The data used by Dobson (1970) consisted of stereophotogrammetr ic 
values of p(x, y) along a line on the wavy surface. Spectra of these values 
were computed from observations made as the wind varied from 7 to 10 knots 
to 1 6 to 1 8 knots. 

The resulting spectra are difficult to interpret. It would be 
necessary to define S(i, m) by means of equations (2.4) through (2.9) over 
five regions defined by circles in the &,m plane. Then it would be necessary 
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to evaluate St**) as given in equation (1.16) for those line sections m 
the upwind -downwind direction by an integration over m broken up into 
pieces corresponding to each spectral range of definition. Since angular 
effects, to be treated later, are not well known, the final result would be 
complicated by its analytical representation of the form used for the angu- 
lar effects. 

The original spectra, despite all of these qualifications, and 
except for some stray constant that cannot be explained, show what would 
be expected on the basis of the proposed wind speed dependent spectra. 
They follow an (I *f 2,5 law instead of an (t *f 3 law. They are wind speed 
dependent in that the spectrum for 16 to 18 knots is appropriately higher 
than the spectrum for 7 to 10 knots. 



S ( k ) [cm 



Fig. 8.1 S(k) versus k for light wind speeds. The graphs are area preserving. 
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Figure 8,2 shows graphs of, and portions of the graphs of S(k) 

versus k plotted on logarithmic scales for u^ equal to 12, 24, 48, 96, 

and 192 cm/sec so that the 19.5 meter wind varies from 3.5 to 34.3 m/s. 

The graphs for u^ = 12, 24, and 48 cm/sec are complete. For u^ = 12, 

S (k) extends to k s 0.35, there is no Kitaigorodskii range, and the 

other three ranges follow. This is the lowest spectrum described by this 

theory, and it corresponds to a significant wave height of 25 cm. For a 

wind slightly lower than this, with no "gusts” above u^ = 12, the capillary 

range would drop by four orders of magnitude. 

For u = 24, S L (k) extends to k “ 0.09, the Kitaigorodskii range 

shows up out to k =0.35, Because of the different behavior of and 

S 4 (k), where one is a function of u^ and the other a function of Dfu^), 

the break in slope of the Leykin-Rosenberg range is pronounced at first, 

then less pronounced, and finally more pronounced. 

The graphs of the spectra for u^ = 96 and u^ =192 are incomplete. 

3 

Only the part of S^(k) for which it is essentially equal to a/2k is shown. 
The Kitaigorodskii range begins at k =0.007 for u^ = 96 and k =0.0018 

for u = 192. A wind of 20.6 m/s (38.5 knots) and waves 9 meters 

* 

(29.5 feet) in significant height can be expected once in a while. 

However, a wind of 34.3 m/s (66.6 knots) and waves 24.9 meters 
(81.6 feet) are not to be expected because the fetch and the duration are 
far too large and never occur. The Sj(k) part of the spectrum would not 
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Fig. 8.2 S(k) versus k with both scales logarithmic for an extreme 
range of wind speeds. 




have the required duration and fetch to be generated. The interval of 
’’small scale isotropic turbulence” covers a large range and could con- 
ceivably be generated in hurricanes. This concept may be closer to reality 
at these high winds than any effort to define a single valued surface, 
q(x,y, t), as shown, for example, in the photographs of waves for winds 
of 60 knots and over in Neumann and Pierson (1966). Just as there is a 
lower bound on u^ at 12 cm/sec for which these results apply, there ls 
an upper bound, less well defined, where the concept of a wind roughened 
sea surface no longer applies. 

Figure 8.3 shows graphs for only a small range of wind speeds. 

They are plotted on a logarithmic k scale and a linear kS(k) scale so 
that area is preserved. Since, from Figure 8.2, S^{k) =30 at k = 0.1, 

kS^ =3 at k = 0.1 and even for an extreme wind a graph of the spectrum 
would hardly move away from the horizontal axis at k = 0.1. Effects of 
higher winds would show on these graphs at k - 0.02, however. 

An important feature of Figure 8.3 is that it shows that the contri- 
bution to the variance of the sea surface elevation comes from wave numbers 
smaller than 0.1 and that for most wind speed wave numbers larger than 
0.01 contribute a very small percentage. 

Slope Spectra 

Figure 8.4 shows area preserving graphs of the slope spectra as 

u varies from 12 to 48 cm/sec corresponding to winds that vary from 
* 

3.5 to 12.3 m/s. Important contributions to the slope variance are pro- 
duced by spectral components with wave numbers that vary all the way 
from 10^ to 30. 
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For the gravity wave -gravity equilibrium range IcS^k) essentially 

rises to the constant value of a/2 for the equilibrium gravity range and 

becomes horizontal until the Kitaigorodskii range rises from it. The con- 
3 3 

’ tributions of the k S^k) and k S^(k) ranges to the slope variance are 
substantial, but the most important contribution comes from the capillary 
range for which equation (8.4) applies. 


o aD(u ) 

k J S 4 (k) = 2^ 


(8.4) 


If U 19 5 were simply a constant multiple of , the spectra would 
not bunch together with increasing u^ on the left and the contribution of 
Sj(k) from k = 0 to k = k^ to the slope variance would be a constant. As 
a result of equations (3.1) and (3.2), however, the contribution of Sj(k) 
to the slope decreases with increas ing wind speed. The values of k^ can be 

located on the line, a/2, as the points where the Kitaigorodskii range begins. 

Other graphs to attempt to illustrate properties of the slope spectra 

were prepared but are not included in this report, A simple linear plot 

is very misleading and difficult to graph. 

Both Figures 8.2 and 8.4 show how the equilibrium gravity range 

depends on wind speed. The upper end of its range is given by k^ and the 

lower end depends on how closely equation (2.5) is required to approach 
3 

a/2k . The lower end of the range also decreases with increasing wind speed. 
Curvature spectra 

5 

Figure 8,5 shows area preserving graphs of k S(k) on a linear scale 
versus k on a logarithmic scale so as to illustrate the range of wave 
numbers that make important contributions to the curvature of the sea sur- 
face. The range is from k = 1 to k = IOtt, corresponding to wavelengths 
of about 6 cm to 0.2 cm. The peak varies from k = 2rr to 4 tt, or wave - 
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Fig. 8.3 kS(k) versus logk for a small range of wind speeds. 
The graphs are area preserving. 
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lengths of 1 cm to 0.5 cm. 

An important feature of this spectrum is that higher moments 

_ 9 

exist. The k behavior for the Cox viscous cutoff range ensures that 
moments up to the seventh can be found. Whether or not they have any 
real physical meaning is, of course, a question. Quantitative measure- 
ments of sea surface curvature do not seem to be available. 

The three figures that show area preserving graphs of the elevation, 
slope, and curvature spectra (Figures 8.3, 8.4, and 8.5) should be com- 
pared. The elevation and curvature spectra each cover about two orders 
of magnitude on the wave number scale, but these two ranges, are sepa- 
rated by more than one order of magnitude. In contrast, the slope spectra 
cover five orders of magnitude with important effects from waves from 
100 meters long to one centimeter long. 

Figure 8.6 shows graphs of curvature spectra (k 4 S(k) versus k) 
on linear vertical and horizontal scales. The only important ranges, as 
far as curvature is concerned, are the capillary and viscous cutoff ranges 
and all the other spectral ranges are to the left of k =0.92 where the 
curves stop. The left hand portion of each spectrum is simply a straight 
line with a slope that depends on D(u^), which, of course, increases with 
increas ing wind speed. 

The curvature spectra given by Figures 8.5 and 8,6 show the 
spectra of the "mean" curvature, which is an invariant property of the 
surface with reference to rotation, as defined by Longuet-Higgins (1962). 
Many other moments of S(£,fi7) up to the fourth can be combined, as shown 
by Longuet-Higgins, to obtain the curvatures in various directions and other 
invariant properties of curvature with reference to rotation. 
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8.6 k 4 S(k) versus k. The graphs are area preserving. 




Freq uency spectra 

Figure 8.7 shows the spectra that result from transforming the 
equations for S^k) through S 5 (k) (equations (2.4) through (2.9)) from 
wave number space to frequency space by means of the general relationship 
given by equation (1.15). Straight lines in wave number space become 
curved lines in frequency space when capillary effects become important. 
The spectra cover 12 orders of magnitude for the range of u^ used and 
3.5 orders of magnitude over frequency. The two small segments near 
w = 5 show the results of Saenger (1972) for light winds of 6 knots 
(3.08 m/s) and 3 knots (1.54 m/s) which correspond to u^ = 10.9 and 
6 cm/sec as discussed in the section on the gravity wave -gravity 
equilibrium range. 

Kondo et al. (1973) have measured frequency spectra for a 
fetch of 40 km and a wind at 10 meters of 6. 1 m/s and for a fetch of 
30 km and a wind of 12. 4 m/s. They report that "Our field results 
do not indicate a capillary saturation limit, at least within the range 
of winds encountered (about 1 - 1 6 m/sec at 10m)". Two spectra are 
shown that cover 9 orders of magnitude in S(OJ) over the range from 
W = 6 to W=200. The similarities to Figure 8.7 are striking except 
that the spectra are higher and grow by more than the results of 
this study indicate they should. Waves in the capillary band seem to 
be too high to be realistic. 
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ANISOTROPY, SUN GLITTER AND RADAR BACKSCATTER 


Introduction 

The data on the form of F(k, $) such that S(k, $) = S(k)F(k, $) 
can be found, now that S(k) has been defined and Illustrated over its 
five wave number ranges, comes from a wide variety of sources. It was 
originally intended to separate radar backscatter results from actual 
wave measurements in discussing F(k, $), but the radar measurements 
provide additional useful information on F(k, $) that can be incorporated 
into an understanding of the angular properties of the spectrum. Com- 
bining various sources that provide information on angular effects pro- 
vides some information about F(k, $) for high wave number ranges, but 
full details are lacking. 

The available data sources and data analyses that provide some 
information on angular effects are those of Cote' et al. (I960), Longuet- 
Higgins, et al. (1963), Toba as given in Figures 5.8, 5.9, and 5.10, Cox 
and Munk (1954) (on sun glitter), Leykin and Rosenberg (1971), Bradley 
(1971) (on the MSC radar backscatter measurements), Valenzuela, 

Laing and Daley (1971) (which provides a partial verification of the 
form of S(k)), and unpublished radar backscatter data provided by 
Dr. J. W. Wright. 

Taken altogether the data provide a fairly consistent description 
of some of the angular effects. They also indicate the need for further 
measurements and for careful interpretation of radar sea return measure- 
ments for winds from 3.5 to 10 m/s, or so. For example, for light winds 
the strongest sea return may come from angles that are 30° or 45° 
away from the upwind downwind direction. 
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Further work on the form of either S($,w) or S(k, <£) is needed. 
This is a difficult problem, but theories and techniques are available, 
such as those of Stilwell (1969). 

Di fferences between measurements under laboratory and natural conditions 

The major difference between laboratory and natural measure- 
ments is that there are atmospheric turbulent scales that are not repro- 
duced in wind-water tunnels. The photographs in Figures 5.8, 5.9 and 5,10 
show for nominal winds of 5.1 m/s that the capillary waves are not travel- 
ing in the wind direction but at angles of 30 to 45° to each side of the wind 
direction. This feature of the wave spectrum was predicted by the reson- 
ant theory of wave generation of Phillips ( 1957). The wind in the wind water 
tunnel blows dominantly down the tunnel with little lateral variation. 

In contrast, winds over the ocean vary from side to side about 
the average wind direction by amounts that would smear any well-defined 
two -peaked form for the angular variation of the spectrum for a given 
wave number, especially if the spot illuminated by the radar is large. 

Some of the effects of turbulence on radar sea return measurements were 
discussed by Pierson and Moore (1972), but this particular effect was not 
mentioned. 

Turbulence is studied in terms of the horizontal, transverse 
and vertical components of the instantaneous wind. The transverse com- 
ponent would have an important effect on the capillary range by producing 
small local patches of short waves traveling in different directions for 
each patch. Studies of angular effects for S(k, $) for both natural con- 
ditions and under laboratory conditions are needed to clarify these features 
of the problem. 
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Theoretical considerations 


For wave numbers corresponding to k > kj equation (2.3) 
needs to be replaced by equation (2.14) which is repeated here as 
equation (9.1) as defined over - ^ ^ and by zero otherwise. 

i r N i 

F(k, u + ) = — 1 + E a n (k, u^) cos 2n$ 

L n= 1 _ 


(9.1) 


This form is inconvenient if the spectrum is required in the 

form S Moreover, so little is known about the a that it will 

n 

only be possible to determine something about a t and a 2 . Equation 
(9.1) is therefore specialized immediately to equation ( 9 . 2 ). 


F(k, u^) = — [l + aj (k, u^) cos 2$ + a^(k, u^) cos 4$] (9.2) 

This can be rewritten as equation (9.3). 

F(k,S; Uj!c ) =|[1 + a 1 (k,u s|| )[2 co S 2 $ - 1] 

+ a 2 (k, u^)^ cos 4 $ - 8 cos 2 $ + 1] ] 

= | [(1 " a 1 (k,u j{c )+a 2 (k,u )!c )) 

+ (2a 1 (k, u^) - 8a 2 (k, u^)) cos 2 $+8a 2 (k, u^) cos 4 $] (9.3) 


Equation (2.3) is a special case of this more' general ex- 

• 4 i 

press ion. For small k, for example, = 3 - and a^ = j . 

Each of the various spectral ranges can now be expressed as 

S(j l,m) where 

7 7 - 

S (uSmV) 2 ,1 lm 

S. (I, m) dro = r— 1 — *F((£ +m ) 2 , tan (^))didw (9.4) 

U + m *)* 1 

because 
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(9.5) 


S.(k)F(k,*) 


k dk d$ 


The results for the Kitaigorodskii range, the Leykln-Rosenbe rg 
range and the capillary range are given by equations (9.6), (9.7), and 
(9.8). 

S 2' 


au 


S ,U,m) = r — — — — 2 7/4 * 


(9.6) 


2k u + m ) 

2 *m ' 


S 3 (A,m) 


S 4 (ji,m) 


aD(u^) 

2k 3 P(j2 + ». 2 ) (4 - p * 77 ' 


aD(u ) 

i—y • F(j l,m) 

2 (l z + m z ) 


F U,ro) 


(9.7) 


(9.8) 


F(i, m) is given by equation (9.9). 


F(t,m) = ^ 


2 2 


(1 - a 1 (( ^ + m Z ) Z , U^) + a 2 ((I + to ) 2 ,%) 


+ (2a 1 ((£ 2 +m 2 f » u *)" 8a 2^ +,R ) 2,u ^ “2 2 

* jb "t TFl 


2 . . 1 A 


+ 8a ? (U +772) 2 ,u)— - — 7 

2 v (Hr + m Z ) 


(9.9) 


For w polarization, and upwind downwind conditions, Valenzuela 
et al. (1971) give the following equation for 0-^(9; u^, k Q ) after correction 
fo r a spectrum defined only over half a plane in i,-OT space. 

0 -° (0;u A ,k ) = 8TTk 4 cos 4 0|G(e,0)| 2 S(2k sine, 0) (9.10) 

w #0 o u 

For crosswlnd conditions, a similar equation results in which 
S(2k Q s in 0, 0) is replaced by S(O,2k o sin0) . 
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In equation (9. 10), | G(e » 0) | 2 is given by equation ( 9 . i 1 ). 


G( e , 0)| 


(e ~ l)[e(l + sin^Q) - sin^O] ,2 
[c cos 0 + Jt - sin 0 


(9.11) 


Equations (9.6), (9.7), (9.8), and (9.9) can be substituted into 
equation (9.10) so as to obtain expressions for cr^ v for each range of 
definition. For the capillary range and upwind -downwind conditions, the 
result is equation ( 9 . 12 ) which is independent of wave number except for 
a j( 2 k Q sin 0 , u^) and a 2 ( 2 k Q sin 0 , ). 


( 0 ;k ,u.) = 
vv o 


aD(uJ 


I G(c, 0 ) 1 2 (cote ) 4 


[! + a 1 ( 2 k o s in 0 , u^) + a 2 ( 2 k Q sin OjU*)] 


( 9 . 12 ) 


For crosswlnd conditions in the capillary range, the result is 
equation (9.13). 

aD(u ) ~ . 

<r i v ( 0; k 0 . 4 — | G(c , 9) | (cot 0 ) • 


[l _ a j ( 2 k Q sin 0 , u^) + a 2 ( 2 k o sin 0 , u^)] 


(9.13) 


Angular effects for the gravity wave-gravity equilibrium range 

Equation (2.3) appears to be a good representation for angular 
effects in the gravity wave -gravity equilibrium range for 0 < k < kj for 
the idealized conditions discussed in the introduction. Ask approaches 
kj , only the second term is important. Near the peak of the spectrum 

4 

the (cos $) term dominates. This equation is a modification of an 
equation given by Cote et al. (I960) to permit some simplifications in 
planned applications to studies in radar backscatter, which will use the 
recent results of Jackson (1972). 


86 


The use of the series of buoys developed at the National Institute 
of Oceanography in Great Britain as described in Neumann and Pierson 
(1966) and by Longuet-Higgins et al. (1963) has resulted in an alternate 
empirical expression for the portion of the spectrum for fully developed 
wind seas of the form 

F(k,3>) = C q [cos($/2)] q(k) (9.14) 

such that the integral from -it to tt is one and the expression is defined 
over the full angular range. 

Near the peak of the spectrum q is typically 20, and it de- 
creases as a function of k (or w) to values of about 5 or 10 toward high 
wave number (or frequency). It is interesting to note that 

(cos $/2) 20 s (cos $) 4 (9.15) 

, , -.,,20 ,-10 

For example, the left hand side at $ = 'rr/2 equals (v2/£) or * 

and the right-hand side is zero. At $ = -rr/3, (cos $/2) 20 equals (\T3/2) , 

or (3/4) 10 or (9/16) 5 = (17.6)" 1 and (cos $) 4 equals (1/2) 4 = (16) 

Instead of a discontinuity at $ = ±ir/2 for the higher wave 
numbers in this range the actual function probably fairs smoothly to 
zero by being smaller than the values given in the equation inside the 
range of definition and small, but non-zero, for some fraction of a 
radian outside of the range of definition. Or in other terms, owing to 
the lateral gustiness of the wind, there may be some spectral components 
traveling opposite to the average wind direction. 

Sun glitter 

The measurements of sun glitter by Cox and Munk (1954) can 
be used to check the reality of the slope spectrum given by Figure 8.4 
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and to learn about anisotropy in the spectrum for the spectral range 
that provides the dominant contribution to the slope. 

The data tabulated by Cox and Munk (1959) provide information 
on the wind speed at two elevations, on wave height, and on the sea surface 
slope variance in the upwind -downwind (cr^ ) and crosswind (o- *) directions. 

It was shown that the slopes had a nearly Gaussian distribution except that 
the upwind -downwind distribution had a peak toward the upwind direction. 
Slight nonlinearities are therefore evident in the slope properties, but to 
a fair degree of approximation k S(k) does describe the slope spectrum. 

The tabulated data in Cox and Munk (1954) give estimates of the 
significant wave height for each measurement. A plot of wave height versus 
wind speed from these data shows practically no correlation between the 
wave height and the wind speed and a very large scatter. For the con- 
ditions of the measurements, the gravity wave part of the spectrum was 
controlled by conditions other than the local wind such as swell from a 
distance and short fetches. 

This wave height variability had little effect on the slope vari- 
ance, although it probably accounts for some of the scatter. Figure 8.4 
shows that the gravity wave range could be quite different from the ideal- 
ized fully developed wind sea condition as far as slope is concerned and 
yet the effect on slope variance would be small. 

If S(k) is the amplitude spectrum, the total slope variance 
integrated over $ is given by equation (9.16). 

<r s 2 = /k 2 S(k)dk (9.16) 

o 

Figures 9.1 and 9.2 compare the slope variance of this spectral 
model with measurements of Cox and Munk (1954) as a function of both 
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Slope Variance <r c C s cr ( , +cr 



» 

Fig. 9.2 Slope variance versus u^ from the theoretical spectra compared 
with measurements by Cox and Munk. 
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u and U 7 7C - (one of the heights at which Cox and Munk measured the 
winds). Integration of the spectrum from k = .003 cm to k = 300 cm 
(essentially k = 0 tooo) yields variances at least a factor of two high. 

Cox and Munk (1954) stated that the resolvability of their system was 
about 40 cm. Integration to this value yields variances a factor of 
two low. These results illustrate that the optics of their system prob- 
ably imposed a high frequency filter on the sea surface reflection. The 
other line in the figures represents a compromise in the high wave 
number cut-off. Figure 9.2 shows the variance computed by integration 
only to = .9419 cm -1 (i. e. , no capillary or Cox viscous cut-off waves 
are involved). This range of integration still does not produce enough 
slope variance. Figure 9.1 shows that the integration to k = 1.2 cm 
(which includes some of the capillaries) produces a reasonable fit to the 
data. 

Cox and Munk (1954) cite other investigations that yielded sea 
surface slopes a factor of two larger than the values that they obtained. 
The theoretical spectra are producing results of the right order of mag- 
nitude for slope, and the two bounds illustrated bracket the Cox and Munk 
data. It is not inconceivable that the system used by Cox and Munk 
imposed a low pass filter on the slopes at a wave number of 1.2 cm 
which corresponds to a wavelength of 5.2 cm. 

The upwind -downwind slope variance is given by equation (9.17). 
In the absence of more detailed knowledge of the variation of a^ and a ^ 
with wave number, they are assumed to be independent of wave number. 
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2 7 

<r = f 

u J 

- 00 


/ l S(l,m)<U 


■n/Zr-OO - 

dw= / /k cos $S(k)F(k,$)dk 

-tt/2Lo 


d$ 


00 -> tt/2 , 

S/k^kJdk / ±[1 
o -tt/2 17 


+ a^cos 2$ + a^cos 4$]cos 



jv 




(9.17) 


Similarly for the crosswind slope variance, 

i_ a il 
2 4 J 


2 2 
or = <r 
c s 


There is no information on in these two variances and 
2 2 2 

°" u + cr^ = or as it should. 

From Cox and Munk (1954), 

a, 2 = 3 X 10' 3 + 1,92 x 10” 3 U 

and 

cr u 2 = 3.16 x 10“ 3 U 


(9.18) 


(9.19) 

(9.20) 


where U is in meters/sec at 41 feet (12.5 meters). It follows 
that equation (9.21) holds, 

2 . a l 

% 3.16U 1 + 2 

7T = 3+ 1.92 U = 17 (9-21) 

>-T 

and therefore 


„ _ 2.48U - 6 

1 ~ 5.08U + 3 ( 9 *22) 

If a.£ were zero, a wind at 12.5 meters of 2.42 m/s would 
produce isotropy. A wind of 3.5 m/s yields an a 1 of 0. 112. A wind of 
20 m/s yields an a^ of 0,42. The angular variation becomes more pro- 
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nounced with increasing wind speed. 

Wind Water Tunnel Measurements of Slope and Curvature 

Wu (1971) has measured the probability density function for 
upwind - downwind slopes, including skewness, and functions related to 
surface curvature distributions in a wind water tunnel. The slope dis- 
tributions were close to normal, but also showed some non-linear 
properties unexplainable by our model* 

The upwind -downwind slopes measured by Wu (1971) exceeded 
the upwind -downwind slopes measured by Cox and Munk by about 
50%, even without the gravity wave part of the spectrum being present. 
From Figure 9*8 and equations 9* 17 and 9* 18 and from Figure 8 of 
Wu, the total slope would be about 0. 06 at — 27 and 0. 1 14 at — 40. 
These values would plot favorably on Figure 9* 1 i n terms of the total 
slope variance predicted by the spectra of this theory. 

Wu (1971) also found that as the free stream wind speed, 

U ■ increased from about 2 to 3 m/ s, the average radius of curvature 
o' 

decreased from 4 to about 0. 8 cm. For U q greater than 3 m/ s, the 
average radius of curvature decreased slowly and linearly with wind 
speed. The curvature spectra in Figure 8.6 have this property. 
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Sea return data, from the Manned Spacecraft Center Program 

A study by Bradley (1971) has analyzed the results of sea return 

measurements for a vertically polarized fan beam scatterometer at 

15.3 GHz. It was not possible to obtain absolute values of c r° (0) but 

w v ' 

<r° (35°) 

1 A 1 W 

0 log 10 o' 0 (10°) ’ 

which is called a normalized scattering cross section, was obtained and 
tabulated for Missions 119 and 156. The results of Bradley (1971) for 
these two missions were also given by Pierson and Moore (1972) and 
interpreted in another context. 

The measurements of interest here were made at a radar 
wavelength of 2.25 cm at an angle off the vertical of 35°. This cor- 
responds to a capillary wavelength of 1.96 cm, or a wave number of 3.2 cm -1 . 
This water wavelength is near the center of the capillary band and the 
radar backscatter should be describable by equations (9.11), (9.12) and ( 9 . 13 ). 

For upwind and crosswind conditions the results of Bradley 
at 35° can be summarized by equations (9.23) and' (9.24) with U in 
knots at 19.5 meters. ! 

10 log 10 [<r“ v (35°)/c^(10«)] u - -34.42 + 14.93 log 1() U (9.23) 

10 log 10 [o-“ v /35»)/ ( r^ v (10o)] c = -35.51 + 12.62 log 1Q U (9.24) 


From equations (9.12) and (9.13) and equations (9.23) and (9.24), 

^u^lVu _ r 35‘ u _ 1 + a i + a 2 _ 10' 3 * 442 U 1 * 493 

<r 35° c / cr 10 o ^35° 1 " a l +a 2 lO- 3 * 551 U 1 * 262 

L c c 

which reduces to 


(9.25) 


94 


(9.26) 



The left-hand side of equation (9-28) represents the values of 

cr° tabulated by Bradley (1971). The right-hand side is the sum of anum- 

w 2 4 

ber of logarithms. Of these, a, |G(e,0)| (cot 0) are known, and 

<r° (10°) is believed to be nearly a constant but is not known for these 
w' 

measurements. Equation (9.28) can be rewritten as equation (9.29). 


<r° (35°) 

10 lo §10 = K + 10 1 ° S 10 D < tt . , < 1 + a l + a 2> 

vv' u 


The values tabulated by Bradley should differ from 

10 l°g 10 Dfu^fl + aj + a 2 ] 


only by some unknown constant. It is tempting to use equations (9.22) 
and (9.27) to solve for both a : and a 2< However, since the equations 
apply to different ranges, and since the true situation is undoubtedly quite 
complicated, it will be assumed that a 2 is zero. 

Table 9.1 shows how this unknown constant was determined and 
how D*(u )[1 + a j] and D*(u + ) were calculated. A conversion from knots 
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tThe cross section values are 


to meters per second is shown in the second column of the table. The 
value of u^ from a neutral atmosphere is shown in the third column. 

The normalized values given by Bradley are shown in the fourth column. 
The antilogs of the four values tabulated for 16 knots were found to be 
0.0237, 0.0316, 0.0204 and 0.0282. These average to the value 0.026 
and D*(u jj< ) from Figure 5.6 or equation (2.12) is 4.85 for this wind 
speed at 19.5 meters and 1 + a^ equals 1.414. Therefore 0.026x = 

4.84 (1.414) and x = 187 (1.414) for the unknown constants (as far as 
these measurements are concerned) in equation (9.29). The addition of 
24.21 db to Bradley's values ought to produce a number equal to 
10 log 1Q D*(u :4t )(l + ap. This value is shown in the fifth column and the 
appropriate antilog is tabulated in the sixth column. 

The next column shows the value of (1 + a^) calculated from 
equation (9.26) with a^ set equal to zero. The next but last column shows 
the value of D^fu^) obtained by dividing D*(u + )(1 + a^) by (1 + a^). The 
last column shows the results of assuming that a^ is a constant equal 
to 0.414. 

Sixteen of these twenty values of D*^) obtained in this way, with 
the correction for variable a^, where the asterisk designated that it was 
calculated from sea return data and not wave data, are plotted on Figure 
5.6 as given in the section on the capillary portion of the spectrum. The 
lowest four values do not fit on the scale used, and the ones for 6 knots, 
or 3.08 m/s, would not be expected to do so. 

The points for 9.3 m/s ought, of course, to straddle the curve 
for D*(u ) because these points were used to determine the unknown 
constant. The interesting feature of this calculation, however, is that 
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these calculated values of D+fu^) also fit the D(uJ curve at 11.3 and 
15.8 to 17 m/s about as well as the values of D(u^) found from measure- 
ments of the capillary spectrum in wind-water tunnels by Toba and 
Sutherland. 

a ^ so plotted against u^ in Figure 5.2. Some values are 
not shown because the points are too crowded together on this scale. The 
points for u^ between 10 and 20 cm/sec seem to be lower than the cluster 
of seven NYU values and D*(u^) obtained from the radar data is lower than 
the values given by the defining equation for D(u jjc ). This discrepancy will 
be explained on the basis of the behavior of a^ to be described later. 

The distance covered by the aircraft in a typical data run for these 

4 

missions was more than 10 km, or 10 meters, compared to typical wind 

water tunnel fetches of 10 meters. If the values of D*(u^) given in 

Table 9.1 were plotted on Figure 5.3, the figure would have to be extended 

at least 1.5 orders of magnitude to the left and the points would form a 

separate isolated pattern paralleling the wind tunnel points over a com- 

2 1 

pletely different range of (u^ /gF) 2 * This strongly suggests that capillary 
waves do not know the difference between a fetch of 10 meters and a fetch 
of 10 kilometers, or more. 

This result lends support to the argument that capillary spectra on 
the open ocean are quite similar to those measured in wind water tunnels 
and provides an important link between laboratory and open ocean measure- 
ments. It is unfortunate that a primary calibration for <r° was not 
possible in the MSC program. However, the agreement, once the data 
are fitted at 16 knots, is striking. The alternate hypothesis that D(u ) 
on the open ocean differs by say, a factor of ten from D(u # ,) in wind water 
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tunnels leads to great logical inconsistencies. 

It should be noted that U irk - is not linearly related to u so that 
the curvature in the graph for D(u^) in Figure 5.6 need not necessarily 
show that 

«r° ~(U 19 , 5 ) 1-49 (9-30) 

does not hold, as found by Bradley (1971). The range from 7, or 8, m/s 

1 49 

to 18 m/s can be fitted to a U power law. The low values for 6 and 
10 knots in Table 9.1 , which may be in part due to going below the critical 
speed of 3.5 m/s, make the data fit this exponential law approximately 
over its full range. 

Spectra deduced from radar sea return 

In this section, an attempt is made to relate wave number spectra 
as deduced by Valenzuela, Laing and Daley (1971) to the preceding re- 
sults. Radar sea return data were used to calculate wave spectra by means 
of the assumption that (the normalized radar cross section) is given 

by Bragg scattering theory. Points were obtained for the wave spectra (in 

4 7 2 

units of cm ) for a range of more than 10 and for a range of 10 in wave 

number by applying the theory to four different radar frequencies and a 

variety of depression angles. The values that were obtained correspond 

to S(£,0) for upwind down wind conditions for values of £ that correspond 

to 2k sin 8 when k is the wave number of one of the four different radars 
o o 

and sin0 is the angle off the vertical, which was varied between 30° and 60°. 

For a set of observations for which the wind (probably at 19.5 meters) 
varied from a maximum of 23 to 24 m/s to a minimum of 2.5 m/s, W (t , 0) 
was calculated from equation (9.31). 
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V [m/s] 



Fig. 9.3 Values of M as a function of V. 
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(9.37) 


(9.38) 


Figure 9*4 shows S(£,0) plotted versus £ on a db scale 
(or equivalently a double log scale) for 5 equal to 3.5 m/s, 12 m/s, 

and 23.5 m/s. The lowest curve represents the spectrum for the mini- 
mum wind speed for which the capillary waves have a substantial ampli- 
tude. To the left of the £ = 0.35 line, the spectrum shows the saturation 
value of the Phillips equilibrium range for the gravity wave range. There 
is no Kitaigorodskii range. 

The middle curve is the curve for a wind of 12 m/s to correspond 
to the data for a wind between 11 and 13 m/s. The Leykin-Rosenberg 
range continues into the Kitaigorodskii range without much change in slope. 
The Kitaigorodskii range will intersect the line with a slope of -4, and 
labeled 1.71 X 10 -3 U' 4 ), at a point just to the left of -12 db. 

The upper curve is for a wind of 23.5 m/s (for the range from 
23 to 24 m/s). The three different slopes for S 2 U), SjU), and S 4 (i) 
are clearly evident. 

The data from Valenzuela et al. (1971), as given in Figure 2 of 
that paper, do not fit these curves. The data appear to be about 11 db 
too low, a factor of about 4 tt, after reduction by 3 db to account for the 
difference in definition. A total discrepancy of about 11 db would yield 
a factor of 12 so that the waves over the range of definition would be only 
29% as high as calculated from wind-water tunnel measurements. 
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Fig. 9.4 S(£, 0) versus 1 in db and on logarithmic scales for winds of 

3.5, 12, and 23.5 m/s compared with data from Valenzuela 
et al. (1971). 





The 11 db discrepancy does not seem to be explainable on the 
basis of any inconsistencies in the basic definitions of either the wave 
spectrum or the normalized radar cross section. Guinard et al. (1971) 
have given details on how the radar data were processed and definition 
of <r° given by their equation 3 is the appropriate one for the calculation. 
Some other definitions do not include the 4 tt. 

The 1 1 db discrepancy appears to be present in the originally re- 
ported values of <r^ (Daley et al. , 1970). After adding +3.2 db to account 
for differences between median and average values and for one-way versus 
two-way beam widths as Valenzuela et al. (1971) did also,* the tabulated 
values at 9 = 60° are about 11 db lower than the values measured by Wright 
in a wind-water tunnel (Table 9.3) for nearly the same radar frequency. 

O (u.^) [ 1 + a^] calculated from equation (9.39) (with 19.42 replacing 20.32) 
is 0.85 for 11 m/s and 1.91 for 23.5 m/s. 

However, once 8 db are added to the values plotted in Figure 2 
of Valenzuela et al. (1971) (3 db are already there), the result is shown by 
the coded points on figure 9.4. The small squares and the large squares 
are for the 23 to 24 m/s winds; the circles and double circles are for the 
11 to 13 m/s winds, and the x's are for the 2.5 m/s winds. Since the M 
values for the 17.5 to 19.5 m/s and 18.5 to 20 m/s cases lie between 
the values for the two cases used, similar results would be expected for 
these measurements also. 

Of the thirteen points for the 23 to 24 m/s winds, nine are quite 
close to the theoretical curve. The six points in the capillary range fit 
an H~ 4 line very well. Three of the remaining points fall close to the 
Kitaigorodskii range. The sharp break in slope from one side of the line 
l = 0.92 to the other is quite evident. 

^Personal communication from Dr. j. W. Wright. 
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Of the fifteen point for the 1 1 to 13 m/s winds, twelve points 

come close to the theoretical curve* The capillary range is clearly one 
_4 

with an & slope. Four points come close to the Kitaigo rodskii range 
near l = 0. 1* 

The large squares and double circles do not come close to the 
theoretical curves. For l =0.3 and 0.5, the poorly fitting points are 
paired and they are the right distance apart so that a shift of 3 db would 
put them on the correct line. 

The x's are for a wind of 2.5 m/s. They fall below the curve 
for 3.5 m/s. They should be expected to be even less than shown. A 
few isolated patches of wind above 3.5 m/s could produce the observed 
results for the points plotted. The x T s in the original paper for low wave 
numbers are up above the curve for a wind of 12 m/s and do not agree 
with any particular theory. 

It appears that the values of the slope found by Valenzuela et al. 

(1971) are fictitious combinations of values composited out of the ranges 

given by equations (9.6), (9.7) and (9.8). Apart from a discrepancy of 

11 db, the data fit the theory given in this paper quite well since 75% 

of the points for winds of 12 and 23.5 m/s are close to the theoretical 

curves. The data also show clearly the -4 slope at the high wave number 

end and the -3.5 slope at the low wave number end. 

In view of these results, the other results given in this study, 

the other conclusions of Valenzuela et al. are very much to be questioned. 

— c / v 

The conclusion that for light winds, <r° is proportional to e , based 
on Pierson and Moskowitz (1964) , is doubtful since the exponential term 
in Pierson and Moskowitz only affects very long waves. The behavior 
with o-° proportional to U V and 2v given by 0.455 for X band (8.9 GHz) 


106 


and 0.619 for C band (4.455 GHz) is also much too low because these 
bands are in the capillary and Leykin-Rosenberg ranges, which are 
strongly wind speed dependent. If Bragg scattering provides a valid sea 
return theory for angles off the vertical of 20° or 30°, or more, than 
one would expect, for radar wave numbers greater than aboutone, that 
the wind speed dependence would be described by Fig. 5.6. 

R a da r backscatter ing data obtained in a wind water tun nel_ 

Dr. J. W. Wright has made some absolute measurements of o- yv 
as a function of wind speed for a value of 0 equal to 60° and a radar 
frequency of 9.375 GHz in a wind-water tunnel at a fetch of 39 feet. The 
results were given in preliminary form, as provided by Dr. Wright in 
Pierson and Moore (1972). They are subject to modifications as the study 
by Dr. Wright continues, but these preliminary results seem to be ex- 
plainable in an absolute quantitative sense by the capillary portion of the 

spectrum proposed by equations (2.8) and (9.8). 

The frequency of 9.375 GHz and the angle of 60° yield the result 
that k = 2k sin 60° = 3.40 and the water wavelength involved is 1.85 cm. 
This wave number is in the capillary range defined by equation (2.8) and 
shown in Figure 8.2. 

For this frequency, e = 54 -i 38.5, | G( c , 9) | 2 = 20.32, and tpot 0) 

at 60° equals 1/9. Since the wave number is fixed, aj and a 2 in 
equation (9.12) can only be functions of u^ . 

Equation (9.12) therefore becomes equation (9.39). 

% V (U ^ = T5 (20 ‘ 32) D(U * )[1 + a l (u * } + a 2 (u * )] (9 ‘ 
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The only remaining problem is to relate the free stream wind 
speeds (converted from feet/min to meters /sec), as measured in the 
wind-water tunnel to u^ values. Values of the free stream wind speed 
(UJ reported by Toba and Sutherland (1967) and measured at N. Y. U. 
for various fetches are given along with the resulting values of u in 
Table 9.2. If it is assumed that similar wind profiles would be observed 
in Wright's facility the values of u^ can be extracted from these data. 

A plot of u^ versus using these values has considerable 
scatter. The extrapolation of the wind profile toward the surface and 
the estimate of its slope to find the stress at the surface is a source of 
this scatter. Quite a bit of the scatter of D(u ) plotted against u in 

V >Jc 

the section on the capillary wave spectrum is probably attributable to this 
scatter. 

To determine the values of u^ to use with Wright's data, u^ was 
plotted against from Table 9.2. A best fit curve and curves for upper 
and lower bounds on the scatter of points were drawn. The values of U 

0O 

from Wright were entered on the graph and a lower bound, best fit, and 
upper bound value for u was determined for each U . 

™ oo 

The result is shown in Table 9.3. The first column is U and the 

oo 

next three columns are the lower bound, best fit and upper bound values 

of u A . The fourth column shows the measured values of <r° obtained 
' w 

by Wright. 

These values are plotted in Figure 9.5. The best fit line is the 
solid line and the two bounds are the dashed lines connecting the x's. The 
break in the curve does not come at u^ =12; it could range from u* = 5 
to u^ =11 with the best fit corresponding to u^ = 8. 
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Table 9.2. Values of (m/s) (the free stream wind speed), 

fetch and u_ (cm/s) from three wind-water tunnels. 


Source 

Fetch (ft) 

U oo 

u * 

NYU 

20.7 

1.46 

5.9 

1 1 

20.7 

2.74 

12.3 

1 ! 

20.7 

4.15 

18.9 

Toba 

32.8 

5.1 

33.5 

11 

44.6 

5.1 

36.7 

1 1 

32.8 

6.1 

39.4 

1 1 

44.6 
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11 
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1 1 

44.6 

7.5 

48.9 

If 
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8.7 
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1 1 

44.6 

8.7 
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1 1 

32.8 

9.7 

64.8 

1 1 

44.6 

9.7 

73.7 

1 f 

32.8 
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81.7 

II 

44.6 
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92.6 

If 

32.8 

1 1.2 

113.0 

1 f 

44.6 

11.2 

102.0 

n 

32.8 

12.1 

112.0 

1 1 

44.6 

12.1 

1 14.0 

Sutherland 28,8 

6.74 

56.1 

n 

43.8 

7.10 

56.7 

1 1 

28.8 

9.64 

80.2 

1 1 

43.8 

10.1 

74.1 

ii 

28.8 

12.1 

102.4 

1 1 

43.8 

12.65 

107.3 

M 

28.8 

15.36 

177.0 

1 ! 

19.1 

15.27 

166.1 
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(The last four columns are In db, ) 



Along with the previously mentioned results of Cox (1958), this 
is an indication that the critical u maybe somewhat smaller than 12. 

T 

The critical value depends on the stability of the air in the tunnel and other 
factors such as the amount of artificial turbulence in the air flow. Perhaps 

the results of previous sections could be extended to u values as low as 

❖ 

9 cm/sec and still give useful results* 

As given in the sixth column of Table 9*3, the circles on Figure 
9.4 show the result of evaluating equation (9.39) for a^ and a^ equal to 
zero using the best estimate of u^. Bragg scattering theory would yield 
these results if the wave number spectrum were isotropic over a half-plane 
in j %-m space. 

The triangles on Figure 9.4 show the results of evaluating equation 
(9.39) for and a ^ = 0. These values are tabulated in the seventh 

column of Table 9.3. 

The isotropic theoretical values fit the data obtained by Wright 
near u >Jc = 25 and 30 cm/sec. From u >1( equal to 30 to equal to 80, the 
anisotropic and isotropic theoretical points bracket the observed values and 
the observed values increase from values near the isotropic values to values 
near the a^ = values over this range. 

Above u^ = 80, the curves diverge. It is understood that there 
were breaking waves in the experiment. Since u^ = 80 corresponds to 
winds of 18 m/s (~36 knots) there may be other important effects over the 
ocean that would modify sea return results. For such high winds, passive 
microwave techniques show promise of being a useful indicator of wind speed. 
It should be noted that the results of Bradley (1972) extend to this range 
also. 
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Between u^ = 12 and u^ = 25, the measured values fall below 
the values computed for an isotropic spectrum. Either the values of Dfu^) 
are too high, or the spectrum is anisotropic with an angular variation such 
that the value of S(f,,0) in the upwind -downwind direction is less than that 
which would be predicted for an isotropic spectrum. 

The second alternative is to be preferred both from the resonance 
theory of wave generation of Phillips (1957) which predicts a bimodal dis- 
tribution of spectral values for light winds and from photographs such as 
Figures 5.8, 5.9, and 5.10 that clearly show dominant trains of capillary 
waves traveling at angles to the wind direction. This same phenomenon 
has been documented many times in the scientific literature. 

The function F(k, u^), defined previously, depends only on u^ 
and $ since k is fixed. It is repeated here as equation (9.40). 

F(k, $ ; u^) = ^ [1 + a^u^) cos 23>+a 2 (u jjc )cos 4$] (9.40) 

Since the values of <r^ v from equation (9.39) are only in the upwind- 
downwind direction, this reduces to equation (9.41). (The t r vanished 
when equation (9*39) was derived.) 

F(0; u^) = i [1 + a 1 ( u !|! ) + a 2 (u :jt )] (9.41) 

With two free parameters to pick up as a function of u^ it is not 
at all difficult to obtain a close fit to the measured values of Wright. For 
example, a ^ and a 2 can be defined by equations (9.42) and (9.43). 
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0 


a 


1 


u - 25 

* 

80 - 25 


3 


12 < u . < 25 

* 

25 ^ u >Jc < 86 


a 2 = 


< 


25 - u 


*Y 


25 - 12 y 0 * 315 12 < u * < 25 


25 ^ u . < 86 
* 


(9.42) 


(9.43) 


These wind speed dependent coefficients were used to evaluate 
equation (9.39) as u^ varied from 12 to 86 and the results are given in 
the last column of Table 9.3. The theoretical values differ from the 
measured absolute values by at most 0.7 db and in general by less than 
0*4 db. 

Figure 9.6 shows the values obtained by Wright connected by 
a solid line, the values for the isotropic assumption as x's connected by 
a dashed line, the values for a^ = ^ and a^ = 0 as x's, and the values 
computed from equations (9.39), (9.42) and (9.43) as circles. The agree- 
ment for the circles is good from u^ = 12 to u^ = 85, corresponding to 
winds at 19.5 meters of 3.5 m/s to about 19 m/s. 

Figure 9.7 shows F(k, u^) for this particular k as a function 
of $ for different values of u^ ranging from 12 to 80. Within an ability 
to resolve such a complicated function by two harmonics, the curves shown 
are quite reasonable as a possible variation of F(k, $;u^). Actually the 
bimodal peak for low u^ values is probably not at plus and minus tt/ 4. It 
is wind speed dependent according to Phillips (1957). Many harmonics 
would be needed to describe F(k, ^/u^) over an adequate range of u^ and 
k, even for the relatively simple conditions found in wind-water tunnels. 
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Measuring <r° (0; u ,k ,$) by varying $ in, say, 10° steps over this 
same range of from zero to 90° could provide some extremely 
valuable results. It is encouraging, however, that reasonable quanti- 
tative agreement between a fairly simple theory for radar backscatter 
and known and reasonably assumed properties of the capillary wave 
spectrum has been found. 

It should also be noted that if a., varies over the ocean in roughly 
the same way that it varies in equation (9.43), the values of Dfu^) would 
be computed to be too low in the analysis of the results of Bradley (1971). 

The triangles in Figure 5.2 for low u^ values would then move up closer 
to the data points obtained by New York University. 

Field measurements by Leykin and Rosenberg 

Leykin and Rosenberg (1971) determined F(k, u^) for one wind 
speed corresponding to a value in excess of u^ = 12. The wind was 5 m/s 
at 20 meters, and the water wavelength was 5 cm. F(2ir/5,$,' 13.4) was 
essentially a constant for - T < $ \ • This is sti11 another indication 

that the spectrum in nature varies from nearly isotropic at low winds to 
a peaked form in the wind direction at high winds. 

Summary of anisotropic properties 

Five different sources have been used to find estimates of the 
anisotropic properties of the spectrum for high wave numbers. Three were 
for particular wave numbers as obtained from the results of Bradley, Wright, 
and Leykin and Rosenberg. The result obtained from Cox and Munk's data 
probably extends down into the Kitaigorodskii range, and the result ob- 
tained from Valenzuela et al. (1971) does extend to that range. 

These five results for a 1 are shown in Figure 9.8 in which 
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Bradley 



Fig. 9.8 aj as a function of u^ from five different sources. 
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equations (9.22), (9.26) (with a., = 0), and (9.37) are all referred to a 

common u and to U 1Q - for a neutral atmosphere, with appropriate 

changes for different anemometer heights and units, A correction for 

atmospheric stability effects has not been made for the Cox and Munk 

result. The dot at U nn c = 5 m/s is from Leykin and Rosenberg’s result. 

iy.D 

The difference between the wave tank result and the open ocean 
results could be explainable on the basis of the lateral component of the tur- 
bulent wind, which would tend to remove the bimodal peakfor low winds. 

For winds of 12 to 20 m/s, or so, a^ is between 1/3 and 1/2. The ratio 
(1 + a ) /(l - a^) is between 2 and 3, and thus crosswind radar backscatter 

values should run 3 db to 4.8 db below upwind -downwind values over this 
range. This, in fact, appears to be what does happen. The next gene- 
ration of radars and radar- radiome te rs should be able to pin down a^ 
to a smaller range. 

For low winds, the different sources yield quite variable pre- 
dictions. This variability is probably a function of atmospheric stability 
and <r /U (see Pierson and Moore, 1972), which can be parameterized 
in terms of atmospheric turbulence theories. 
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CONCLUSION 


The wave number and frequency spectra of a wind roughened 
sea have been described and documented. More and better data are 
needed to improve on these results, in particular with reference to 
the values of and k^ . Also other data not available to us should 
be used to check these various ranges in greater detail, especially 
with reference to the Kitaigorodskii and Leykin -Rosenberg ranges. 

Ways to study the shorter waves in wave number space should be 
developed for open ocean conditions. 

The results of this study can be used as a framework for 
deeper theoretical investigations of the generation and properties of 
the shorter waves. 

Radar backscattering data yield results that agree with Bragg 
scattering theory applied to the wave number spectra proposed for 
the high wave number ranges for laboratory backscatter measurements. 

An unequivocal agreement between radar backscatter measure- 
ments over the ocean and the spectra given in this paper has not been 
established. If the tabulated results of Daley et al. (1970) are correct, 
S 2 (k), S 3 (k), and S 4 (k) would have to be a factor of about 12 (or 4ir) lower 
than the values given. Such a reduction would be difficult to reconcile with 
the substantial overall internal consistency of the proposed spectra and 
the other data sources that were used. In contrast, the results of Kondo 
et al. (1973) are higher over the capillary range and provide support 
for the argument that the capillary spectrum does not decrease with in- 
creasing fetch. 
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